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ABSTRACT 

This  thesis  investigates  various  techniques  for  fore- 
casting the  volume  of  containerizable  cargo  that  flows  into 
the  container  stuffing  station  at  the  Military  Ocean 
Terminal,  Bay  Area,  Oakland,  California.   Cargo  input  data 
is  analyzed  in  terms  of  weekly  cargo  volume  inputs  for  a 
selected  number  of  major  ports  of  debarkation.   The  time- 
series  data  for  these  ports  is  first  tested  for  serial 
correlation.   Based  on  the  affirmative  results  of  the 
serial  correlation  test,  the  following  forecasting  methods 
are  investigated:   the  moving  average,  the  exponentially 
weighted  average,  the  exponentially  weighted  average  with 
trend  adjustment  and  the  exponentially  weighted  average  with 
an  adaptive  response  rate.   By  means  of  statistical  testing 
procedures,  the  "best"  forecasting  method  is  determined. 
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I.   INTRODUCTION 

In  light  of  the  dispersal  of  United  States  forces 
throughout  the  world,  the  shipment  of  the  vast  quantities 
of  needed  supplies  and  material  to  overseas  stations  is  a 
predominant  logistics  problem  for  the  Department  of  Defense. 
Since  the  late  1960's  the  need  for  efficient  and  economical 
handling  of  this  cargo  has  led  the  DoD  to  use  containerized 
shipments  on  commercial  sea-going  carriers  as  the  primary 
routine  transportation  source.   The  retention  of  break  bulk 
type  shipping  operations  has  been  necessary  due  to  the  size 
limitations  of  containers.   Also,  the  fulfillment  of  urgent 
requirements  where  time  is  critical  have  been  met  using  air 
transport.   However,  ocean  containerization  offers  the  most 
cost-effective  method  of  handling  routine  shipments  between 
ports  equipped  to  handle  this  type  of  traffic,  as  are  most 
Pacific  Ocean  ports. 

The  ability  to  fill  or  stuff  a  container  at  the  con-  • 
signor's  warehouse  and  unpack  it  at  the  consignee's  ware- 
house, with  no  intermediate  handling  of  the  cargo,  is  a 
tremendous  asset.   However,  many  shipments  are  not  large 
enough  to  economically  warrant  their  own  vans;  therefore, 
stuffing  with  other  shipments  bound  for  the  same  destination 
may  be  accomplished  at  a  central  location.   One  such  stuff- 
ing facility  is  the  container  freight  division  (CFD)  of  the 
Military  Ocean  Terminal  Bay  Area  (MOTBA)  located  in  Oakland, 
California.   This  facility  stuffs  vans  with  containerizable 

12 


cargo  enroute  to  ports,  mainly  in  the  Pacific  and  the  Far 
East,  and  arranges  transportation  of  the  containers  to  the 
ports  via  commercial  carrier. 

All  commercial  carriers  require  that  bookings  for  vans 
and  space  aboard  vessels  be  made  in  advance  of  the  ship's 
loading  date.   As  of  the  fall  of  1974,  this  lead  time  was 
at  least  three  weeks.   To  effectively  make  these  bookings 
for  an  advanced  date,  a  prediction  of  the  cargo  on  hand  on 
the  date  of  arrival  of  the  vans  is  required.   Normally, 
vans  arrive  at  the  stuffing  station  four  to  six  days  prior 
to  the  loading  date.   Thus,  the  request  for  space  must  take 
into  account  (1)  the  present  volume  of  cargo  in  the  ware- 
house, (2)  the  forecasted  arrivals  of  cargo  during  the 
period  between  the  booking  and  the  van  arrivals  for  stuff- 
ing and  (3)  any  previous  bookings  that  will  deplete  the 
cargo  volume  destined  for  the  port  during  the  forecasting 
period. 

The  general  objective  of  this  thesis  is  to  examine 
practical  forecasting  methods  which  can  be  utilized  to  pre- 
dict cargo  arrivals  at  the  container  stuffing  station.   To 
predict  arrivals  there  are  two  general  approaches  that 
could  be  used.   The  first  approach  is  based  on  notification 
of  shipment  departures  by  consignors,  which  would  provide 
advanced  notice  of  arrivals  at  the  stuffing  station.   The 
other  approach  is  based  on  historical  data,  from  which  pro- 
jections for  future  arrivals  may  be  drawn. 

An  advance  notification  system  does  not  presently  exist 
for  the  majority  of  cargo  arriving  at  the  station;  however, 
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it  could  be  produced  by  requiring  consignors  to  provide  the 
CFD  with  three-week  advance  notification  of  cargo  shipments, 
perhaps  by  post  card.   Use  of  this  approach  would  require 
analysis  of  the  delay  between  the  arrival  of  notifications 
and  shipments  from  each  consignor.   Also,  the  statistical 
characteristics  and  stability  of  the  transit  times  for  ship- 
ments and  notifications  must  be  determined.   It  is  doubtful 
that  most  consignors  have  enough  advance  knowledge  of  a 
shipment's  departure  to  provide  a  timely  and  accurate  noti- 
fication, especially  if  the  departure  is  two  or  three  weeks 
in  the  future.   The  effect  of  the  possible  lack  of  consignor 
cooperation  must  also  be  assessed.   Taking  into  account  the 
resources  that  would  be  required  to  effectively  implement 
an  advance  notification  approach,  it  was  decided  that  fore- 
casting methods  which  rely  only  on  historical  data  for 
future  projections  would  be  more  appropriate  for  a  thesis 
topic. 

The  historical  data  used  here  consisted  of  the  cargo 
arrival  data  for  the  CFD  for  fiscal  year  1974;  it  repre- 
sented the  most  recent  peacetime  data  available.   The  data 
provided  the  volume,  weight  and  date  of  arrival  of  each 
shipment.   Of  the  three  alternative  units  for  forecasting 
(number  of  shipments,  volume  and  weight),  the  volume  of  the 
cargo  was  selected  because  it  is  currently  used  by  the  mari- 
time community  as  the  measure  for  booking  containers  on 
vessels . 

Two  alternative  aggregation  levels  could  have  been  used 
for  forecasting.   One  was  based  on  the  port  of  debarkation 
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(POD) ,  and  the  other  based  on  consignee  destination  of  the 
cargo.   There  are  more  than  1,000  active  consignees  served 
by  the  CFD,  which  would  require  a  very  large  number  of  in- 
dividual forecasts  if  the  consignee  base  were  used.   PODs 
are  the  currently  used  aggregation  levels  for  forecasts  by 
CFD  personnel  because  van  space  is,  by  convention,  booked 
aboard  ships  destined  for  PODs,  not  consignees.   Therefore, 
PODs  were  selected  as  the  aggregation  level.   Analysis  of 
the  data  revealed  16  major  PODs  which  received  over  901  of 
the  shipments  passing  through  the  CFD. 

The  effectiveness  criteria  used  for  the  selection  of  an 
optimal  method  are:   1)  computational  feasibility,  and  2) 
the  ability  of  the  method  to  maximize  forecasting  accuracy 
by  minimizing  forecasting  error.   Computational  feasibility 
is  limited  here  to  those  methods  which  do  not  require  the 
use  of  a  computer  to  routinely  generate  forecasts.   Mean 
absolute  deviation  and  mean  deviation  squared  are  used  as 
measures  of  forecast  errors. 

Forecast  errors  were  computed  over  time  horizons  of 
one,  two  and  three  weeks,  corresponding  to  possible  lead 
times  in  the  booking  process.   In  other  words,  for  a  two- 
week  time  horizon  a  single  forecast  was  used  for  two  suc- 
cessive weeks.   The  deviation  for  the  forecast  was  computed 
by  taking  the  average  difference  between  the  forecast  and 
the  actual  input  volume  for  each  of  the  weeks  covered  by 
the  forecast. 

Generally,  the  selection  among  various  methods  of  fore- 
casting depends  upon  the  statistical  characteristics  of  the 

15 


data.   These  characteristics  are  differentiated  by  the  degree 
of  statistical  dependence  between  successive  observations. 
If  there  is  no  statistical  dependence  between  successive 
observations,  then  ordinary  trending  and  indexing  techniques 
may  be  applied.   On  the  other  hand,  the  presence  of  statis- 
tical dependence  among  successive  observations,  (called 
serial  correlation) ,  implies  a  need  for  methods  which  place 
more  emphasis  on  recent  past  observations  when  predicting 
future  volumes. 

Tests  for  serial  correlation  of  the  data  for  the  major 
PODs  were  made  using  the  Durbin-Watson  d-statistic.   The 
aggregated  results  of  these  tests  indicated  that  the  pre- 
sence of  serial  correlation  was  statistically  significant. 

Figure  1  is  a  graph  of  the  weekly  input  volumes  of  POD 
UL7  for  a  20-week  period  commencing  with  the  nineteenth 
week  of  FY  '74.   This  was  selected  as  a  typical  representa- 
tion of  the  data  for  the  major  PODs.   This  particular  POD 
exhibited  strong  positive  serial  correlation. 

With  the  existence  of  serial  correlation  in  the  data 
confirmed,  the  moving  average  technique  and  three  exponen- 
tially weighted  average  techniques  were  examined  for  pre- 
diction accuracy.   These  methods  were  chosen  for  their 
relative  ease  of  computation  and  ability  to  adapt  to  changes 
in  volume  inputs  over  time. 

Simple  moving  averages  [Ref.  5]  of  from  one  to  ten  week 
durations  were  the  first  methods  explored.  In  each  method, 
the  most  recent  weeks  for  the  designated  duration  are  aver- 
aged and  the  result  applied  to  the  forecast  through  the 
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designated  horizon.   The  three-week  moving  average  was  the 
technique  in  actual  use  at  the  CFD.   The  resulting  data  in- 
dicated that  for  all  time  horizons,  moving  averages  greater 
than  three  weeks  generally  provided  lower  deviations  and 
thus,  better  forecasts  than  the  three-week  moving  average. 
Note  that  the  longer  the  duration  of  the  average,  the  less 
responsive  the  predictions  are  to  recent  changes.   The  data 
indicated  that  over-sensitivity  was  as  undesirable  as  no 
sensitivity  in  most  cases. 

The  next  method  examined  was  the  simple  exponentially 
weighted  average,  [Ref .  5] ,  or  exponential  smoothing  method. 
This  method  used  a  weighting  or  smoothing  constant,  alpha, 
a  value  usually  between  zero  and  one,  to  regulate  the  em- 
phasis on  the  most  recent  observation.   Large  values  of  alpha 
place  increasing  emphasis  on  the  most  recent  observations 
and  less  on  preceding  observations.   A  forecasted  average 
for  an  observation  period  was  computed  as  the  sum  of  the 
observations  in  that  period  times  alpha,  plus  the  forecasted 
average  for  the  previous  period  times  one  minus  alpha.   In 
this  way,  the  new  observation  was  added  to  the  accumulated 
weighted  average,  a  continually  changing  number.   This  fore- 
cast average  was  then  used  as  the  forecast  and  applied  over 
the  next  weeks  of  the  designated  horizon  and  the  deviation 
statistics  accumulated. 

To  discover  the  levels  of  alpha  which  minimized  error 
measurement,  the  levels  of  alphas  were  varied  by  0.01  from 
-0.1  to  1.0,  and  applied  to  the  data  for  each  POD  and  hori- 
zon.  The  alpha  levels  which  produced  minimal  deviation 
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statistics  ranged  from  -0.02  to  +0.77  depending  on  the  POD 
and  time  horizon  specified.   The  optimal  single  alpha  level 
for  all  PODs  was  about  0.35  and  varied  slightly  depending 
on  the  time  horizon. 

The  next  forecasting  method  explored  was  the  exponen- 
tially weighted  average  with  adjustment  for  trends  [Ref .  5] . 
This  procedure  was  examined  to  determine  if  trend  adjust- 
ments would  reduce  forecast  error.   This  procedure  started 
with  the  same  basic  formula  as  the  simple  exponentially 
weighted  average,  however  a  trend  factor  is  also  computed. 
This  trend  is  accumulated  and  smoothed  in  the  same  manner 
as  was  applied  to  the  forecast  averages,  resulting  in  an 
average  trend  adjustment.   The  forecast  volume  to  be  used 
in  the  predictions  was  created  by  combining  the  forecasted 
average  with  the  average  trend  adjustment. 

To  obtain  minimal  errors,  alpha  levels  were  again  varied 
by  0.01  from  -0.10  to  1.00  and  the  method  applied  to  the 
data  for  all  PODs.   The  results  revealed  minimized  devia- 
tions for  alpha  levels  between  -0.01  and  0.38,  again,  de- 
pending on  the  POD  and  time  horizon  specified.   When 
accounting  for  trends,  the  optimal  single  alpha  level  for 
all  PODs  was  0.16,  varying  slightly  for  the  different  time 
horizons . 

The  final  forecasting  procedure  used  was  the  exponen- 
tially weighted  average  using  an  adaptive  response  rate 
technique  to  modify  the  alpha  level  used  from  period  to 
period  to  account  for  fluctuations  in  the  observations 
[Ref.  18].   The  absolute  and  actual  errors,  or  deviations 
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and  actual  errors  between  the  observed  and  forecasted 
volumes  were  separately  accumulated  and  smoothed  in  the 
same  fashion  as  the  forecast  average  using  a  smoothing  con- 
stant, gamma.   A  tracking  signal  was  then  computed  by  di- 
viding the  smoothed  absolute  error  by  the  smoothed  error. 
Deviation  statistics  using  the  simple  exponentially  weighted 
average  were  computed  with  the  alpha  level  for  each  forecast 
set  at  the  absolute  value  of  the  tracking  signal. 

To  evaluate  the  adaptive  response  rate  method,  it  was 
applied  to  the  data  for  each  of  the  PODs  while  varying  the 
time  horizon  and  the  level  of  the  smoothing  constant,  gamma, 
by  0.01  from  -1.00  to  1.00.   Here  the  use  of  negative  con- 
stants indicates  that  the  new  smoothed  error  term  is  created 
by  enlarging  the  old  smoothed  error  and  adding  or  subtract- 
ing a  portion  of  the  new  error  to  it.   The  minimum  deviation 
statistics  for  the  individual  POD  and  various  time  horizons 
were  realized  with  gamma  levels  ranging  from  -1.00  to  0.49. 
The  resultant  optimal  single  gamma  level  for  all  POD  was 
about  -0.05  depending  on  the  horizon  used. 

In  summary,  the  experimentation  with  the  above  forecast- 
ing methods  had  created  deviation  statistics  for  the  follow- 
ing forecasting  procedures: 

1)  Three-week  moving  averages  (maintenance  of  status 
quo)  ; 

2)  Moving  averages  with  a  single,  aggregately  optimal 
duration ; 

3)  Moving  averages  with  the  optimal  duration  for  each 
POD; 
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4)  Simple  exponentially  weighted  averages  using  a 
single,  aggregately  optimal  alpha  level; 

5)  Simple  exponentially  weighted  averages  using  in- 
dividually optimal  alpha  levels  for  each  POD; 

6)  Exponentially  weighted  averages  with  trend  adjust- 
ment using  a  single,  aggregately  optimal  alpha  level; 

7)  Exponentially  weighted  averages  with  trend  adjust- 
ment using  individually  optimal  alpha  levels  for  each  POD; 

8)  Exponentially  weighted  averages  with  adaptive 
response  rates  using  a  single,  aggregately  optimal  gamma 
level ; 

9)  Exponentially  weighted  averages  with  adaptive  re- 
sponse rates  using  individually  optimal  gamma  levels  for 
each  POD. 

Using  these  deviation  statistics,  the  selection  of  the  "best" 
forecasting  method  was  made  subject  to  the  following 
considerations: 

1)  For  operational  consistency,  it  was  assumed  that 
only  one  of  the  methods  would  be  chosen  for  use  with  all 
PODs .   No  attempt  was  made  to  evaluate  the  problem  by 
choosing  the  optimal  forecasting  method  for  each  POD. 

2)  The  deviation  statistics  for  the  one  week  time 
horizon  were  used  in  the  evaluation  as  the  relative  effec- 
tiveness of  the  various  methods  for  each  POD. 

3)  The  mean  absolute  deviation  statistic  was  selected 
as  the  principal  goodness  criterion  since  no  extra  penal- 
ties for  large  individual  forecast  errors  were  intended. 
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In  the  end,  statistical  significance  tests  are  used  to 
reveal  whether  the  meau  absolute  deviations  using  one  fore- 
casting method  were  significantly  less  than  the  deviations 
of  another  method.   The  tests  were  conducted  at  a  signifi- 
cance level  of  0.05.   The  results  indicated  that  the  simple 
exponentially  weighted  average  and  adaptive  response  rate 
methods  were  not  significantly  different  from  each  other. 
However,  both  methods  were  better  than  any  of  the  others. 

Based  on  the  above  results  the  use  of  the  simple  expo- 
nentially weighted  average  was  chosen.   The  selection  was 
founded  on  the  ease  of  computation  of  the  forecast,  as  the 
adaptive  response  rate  method  required  many  more  steps  for 
each  calculation  and  offered  no  better  results.   Addition- 
ally, the  adaptive  response  rate  method  required  a  broad 
range  of  gamma  levels  and  was  characterized  by  multi-modal 
relationships  between  the  gamma  level  and  MAD  for  most  cases 
of  PODs  and  time  horizons.   This  indicated  great  sensitivity 
of  accuracy  to  the  gamma  levels  for  particular  data  sets. 

When  compared  to  the  three-week  moving  average,  the  use 
of  this  simple  exponentially  weighted  average  technique 
provided  4.9%,  8.01  and  10.0%  reductions  in  mean  absolute 
deviations  for  the  one-,  two-  and  three-week  time  horizons 
respectively.   In  addition  to  increasing  the  forecast  ac- 
curacy the  exponentially  weighted  average  is  also  as  easy 
computationally  as  the  three-week  moving  average.   Opera- 
tional benefits  of  the  increased  forecast  accuracy  are  a 
reduction  in  the  average  age  of  cargo  and  volume  of  cargo  in 
the  warehouse,  and  increased  cube  utilization. 
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II.   BACKGROUND 

A.   THE  BEGINNING  OF  CONTAINERIZATION 

While  containers  had  been  used  for  transporting  cargo 
aboard  vessels  as  early  as  the  1820's,  the  modern  era  of 
containerization  began  in  1966  when  Malcolm  McLean  of  Sea- 
Land  Service,  Inc.,  began  loading  containers  aboard  ships 
bound  for  Europe.   The  economics  of  containerization  soon 
became  apparent.   By  the  end  of  1968,  one-fourth  to  one- 
third  of  shipments  on  major  U.S.  routes  were  containerized. 
American-flag  operators  as  well  as  their  foreign-flag  com- 
petitors began  to  invest  heavily  in  the  containerized  sys- 
tem.  This  system  includes  not  only  containers,  but  also 
specialized  ships,  container  terminals  and  road  equipment. 

Within  the  industry  the  container  came  to  be  called  the 
greatest  thing  in  packaging  since  the  paper  bag.  Most  con- 
tainers are  simple  boxes  of  plywood  or  aluminum  and  are 
manufactured  in  standard  sizes:  approximately  8  feet  wide, 
8  feet  deep  and  various  lengths  between  20  to  40  feet  long. 
The  container  volume  capacity  ranges  from  1,100  cubic  feet 
to  2,400  cubic  feet,  and  weight  capacity  ranges  from  38,000 
pounds  to  46,000  pounds  depending  on  the  size  of  the  van. 

Containerization  required  the  shipping  industry  to 
adopt  an  entirely  new  line  of  thinking  and  provided  the  in- 
dustry with  many  phases  of  economic  improvement:   more  ef- 
ficiency in  cargo  and  terminal  handling,  a  reduction  in 
merchandise  damage  due  to  weather  and  handling,  a  reduction 
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in  terminal  facility  requirements,  a  savings  in  reduced 
special  packaging  for  export,  a  reduction  in  pilferage  and 
other  losses,  the  ability  of  manufacturers  to  operate  with 
reduced  inventory,  and  the  benefit  of  computerized  inven- 
tory-control systems. 

Containerization  has  cut  shiploading  and  unloading  time 
from  as  much  as  eight  days  to  as  little  as  twelve  hours. 
The  varied  forms  of  containerization  are  enabling  cargo 
ships  to  spend  as  many  as  nine  out  of  every  ten  days'  pro- 
ductivity at  sea  instead  of  half  their  time  in  port  -- 
causing  some  shipping  lines  to  have  alternating  crews  so 
that  the  men  can  spend  some  time  with  their  families.   A 
typical  containership  is  loaded  with  some  700  vans,  seven 
deep  and  six  across,  and  is  on  its  way  at  22  knots. 

The  new  era  has  resulted  in  transmodal  systems  to  move 
cargo  via  land,  sea  and  air  with  computerized  efficiency. 
In  some  cases,  overall  transport  costs  are  down  30-40%  as 
compared  to  conventional  break  bulk  cargo  handling. 

Sea-Land  Service,  Inc.,  has  become  the  biggest  container 
operator  in  the  world  and  is  one  of  the  two  largest  U.S. 
steamship  lines.   Sea-Land  has  five  33-knot  containerships 
with  capacity  for  1,082  units.   With  other  carriers  taking 
advantage  of  container  economics ,  the  United  States  flag 
carriers  have  become  the  world's  largest  containership 
fleet. 

The  port  of  Oakland,  California  is  now  the  leading  West 
Coast  cargo  shipping  terminal  --  mainly  due  to  containeriza- 
tion.  The  port  handles  45%  of  all  container  shipping  on 
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the  West  Coast  and  is  the  second  largest  container  port  in 
the  world  --  New  York  being  the  largest. 

According  to  James  L.  Reynolds,  President  of  the 
American  Institute  of  Merchant  Shipping  (AIMS) ,  by  1976 
cargo  moving  to  and  from  U.S.  ports  is  expected  to  rise  to 
270  million  measurement  tons  (a  unit  of  volume  representing 
40  cubic  feet)  plus  7.5  million  in  military  cargo  and  an- 
other 20.5  million  in  government  impelled  cargo  such  as 
foreign  aid  shipments. 

Containerization  has  been  recognized  as  a  system  that 
will  best  serve  the  needs  of  the  public  and  the  transporta- 
tion industry.   The  container  revolution  stems  from  the  im- 
mense transportation  capability  gained  from  moving  goods 
swiftly,  safely  and  reliably. 

B.   CONTAINERIZATION  AND  THE  DEPARTMENT  OF  DEFENSE 

In  recent  years  the  United  States  Department  of  Defense 
has  used  commercial  containerization  as  the  primary  means 
of  transporting  its  general  ocean  cargo.   Substantial  cost 
savings  from  the  use  of  containerships  in  lieu  of  using  the 
traditional  break  bulk  cargo  ships  have  been  recognized. 
Yet,  there  still  remain  problems  for  transportation  offi- 
cials within  DoD  to  reckon  with. 

Figure  2  illustrates  the  portion  of  the  transportation 
system  within  the  Department  of  Defense  to  which  this  thesis 
addresses  itself. 

Orders  are  placed  by  DoD  overseas  commands  known  as 
consignees.   Various  vendors  and  military  depots  known  as 
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consignors,  fill  these  purchase  orders  which  the  consignees 
need  to  execute  their  missions.   It  now  becomes  the  job  of 
the  DoD  transportation  managers  to  deliver  these  items  to 
the  consignee. 

Cargo  leaving  the  consignor  for  the  consignee  may  be 
transported  across  the  ocean  in  either  container  or  break 
bulk  ships.   Conventional  ships  transport  break  bulk  cargo 
to  the  requesting  overseas  command.   If  the  goods  are 
ordered  in  large  enough  quantities  and  are  containerizable , 
the  consignor  may  place  the  cargo  in  containers  at  his  ware- 
house.  This  is  known  as  source  stuffing.   The  sealed  con- 
tainer is  then  shipped  to  the  POE  where  it  has  been  booked 
aboard  a  containership  for  overseas  transport  to  the  con- 
signee.  This  type  of  operation  provides  the  most  rapid  and 
economical  form  of  ocean  shipping. 

Another  situation  occurs  when  container  stuffing  sta- 
tions, at  the  POE  or  inland,  receive  quantities  of  contain- 
erizable cargo  for  consolidated  stuffing  into  containers. 
At  the  stuffing  station,  cargo  for  particular  consignees  is 
stuffed  only  after  meeting  volume  and  weight  restrictions. 
A  commercial  shipping  company  then  receives  the  containers 
from  the  stuffing  station  for  lift  aboard  ccontainerships 
and  overseas  transport. 

At  the  Port  of  Debrakation  (POD)  single  consignee  con- 
tainers are  sent  directly  to  the  respective  consignees. 
Containers  with  cargo  for  more  than  one  consignee  are  trans- 
ported to  designated  break  bulk  stations  (BBS) ,  which  may 

itself  by  a  consignee,  where  they  are  broken  down  for  further 
shipment  to  individual  consignees. 
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C.  WESTERN  AREA,  MILITARY  TRAFFIC  MANAGEMENT  COMMAND 
Western  Area,  Military  Traffic  Management  Command  (WAMTMC) 

is  a  jointly-staffed  field  organization  under  Headquarters, 
Military  Traffic  Management  (MTMC) ,  Washington,  D.  C.   WAMTMC 
is  responsible  for  transporting  export  cargo  of  the  Depart- 
ment of  Defense  and  other  government  agencies  to  overseas 
commands.   It  has  been  located  at  the  Oakland  Army  Base, 
Oakland,  California,  since  being  established  in  February 
1965.   MTMC  is  a  major  field  command  of  the  Department  of 
the  Army. 

WAMTMC  is  responsible  for  transportation  management  of 
domestic  and  export  shipments  in  14  western  states:   Idaho, 
Arizona,  California,  Montana,  Nevada,  Oregon,  Utah,  New 
Mexico,  Washington,  Colorado,  Nebraska,  South  Dakota,  North 
Dakota  and  Wyoming.   One  of  the  three  Military  Ocean  Termi- 
nals operated  by  WAMTMC  on  the  West  Coast  will  receive  ex- 
port cargo  from  these  14  states  which  is  destined  for 
shipment  to  overseas  installations.   The  three  terminals 
operated  by  WAMTMC  are:   Military  Terminal  Unit,  Pacific 
Northwest,  Seattle,  Washington  (PNW) ;  Military  Ocean  Termi- 
nal, Bay  Area,  Oakland,  California  (MOTBA) ;  Southern 
California  Outport,  Long  Beach,  California  (SCO).   During 
peacetime,  WAMTMC  also  directs  the  flow  of  all  air  cargo 
in  the  Military  Airlift  Command's  transportation  system 
through  the  Military  Airlift  Clearance  Authority  (MACA) . 

D.  MILITARY  OCEAN  TERMINAL,  BAY  AREA 

Military  Ocean  Terminal,  Bay  Area  is  the  consolidation 
of  the  Army  and  Navy  terminal  facilities  in  the  San  Francisco 
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Bay  Area.   It  represents  the  largest  operating  unit  of 
WAMTMC  with  terminal  facilities  at  the  Oakland  Army  Base 
and  the  Alameda  Reefer  Facility.   Seven  deep  water  berths 
are  under  its  control  as  well  as  an  84-acre  tidewater  con- 
tainer stuffing  area  located  at  the  Oakland  Army  Base. 
Also  under  MOTBA's  control  are  the  ports  of  Stockton, 
Sacramento,  and  Eureka. 

E.   TIDEWATER  CONTAINER  STUFFING  STATION 

The  Container  Stuffing  Station  (CSS)  operates  under  the 
Container  Freight  Division  (CFD)  at  MOTBA.   Containerizable 
cargo  is  loaded  into  the  shipping  containers  at  both  instal- 
lations.  These  containers  are  then  transferred  to  the  com- 
mercial shipping  companies  for  shipment  to  overseas  commands. 

The  CSS  is  operated  under  civilian  contract.   Container- 
izable cargo  is  placed  in  warehouses  until  containers  become 
available  for  stuffing.   In  most  cases  commercial  containers 
are  used  for  stuffing,  although  military  vans  are  sometimes 
used.   The  contractor  has  certain  restrictions  on  container 
stuffing.   Some  examples  are  cargo  mix,  container  utiliza- 
tion, and  consignee  mixing  limitations.   In  addition,  volume, 
weight  and  size  of  containers  vary  among  the  various  com- 
mercial carriers.   There  are  16  types  of  containers  presently 
in  use;  and,  in  general,  each  company  can  accommodate  only 
those  containers  which  have  been  specifically  designed  for 
its  vessels. 
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F.  MILITARY  SEALIFT  COMMAND,  PACIFIC 

The  Military  Sealift  Command  >(MSC)  is  an  agency  that 
coordinates  the  various  ocean  transportation  activities  of 
the  Department  of  Defense.   Within  the  continental  United 
States  (CONUS)  the  relationship  between  the  Military  Sealift 
Command,  Pacific  (MSCPAC)  and  WAMTMC  is  very  close.   In 
order  to  get  its  cargo  aboard  break  bulk  or  container  ves- 
sels, WAMTMC  notifies  MSCPAC  of  the  amount  of  space  required, 
the  number  of  containers  needed  by  the  shipper  and  the  size 
of  the  containers.   MSCPAC  books  the  space  aboard  a  vessel 
and  notifies  WAMTMC  of  the  vessel  sailing  date,  lift  date 
for  the  containers  and  the  number  of  containers  booked. 
This  procedure  is  carried  out  whether  the  shipper  is  a  ven- 
dor or  WAMTMC  is  the  shipper. 

G.  MTMC  MANAGEMENT  INFORMATION  SYSTEMS 

MTMC  employs  two  management  information  systems  in  order 
to  control  the  flow  of  cargo  within  CONUS  that  will  be 
eventually  shipped  to  overseas  commands.   Both  systems 
utilize  a  Burroughs  5500  computer.   The  two  systems  are  the 
Mechanized  Export  Traffic  System  (METS)  and  the  Surface  Ex- 
port Cargo  System  (SURS) . 

METS  monitors  the  DoD  cargo  flow  from  its  origination 
at  the  consignor's  warehouse  until  its  arrival  at  a  state- 
side Port  of  Embarkation,  i.e.,  one  of  WAMTMC ' s  three  ocean 
terminals.   The  consignor,  which  may  be  a  commercial  vendor 
or  a  military  depot,  thus  provides  the  basic  input  into  the 
system.   This  input  includes  shipment  data,  commodity 
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information,  transportation  data,  destination,  etc.   Im- 
properly formated  documentation  is  frequently  received  from 
the  commercial  vendors  requiring  the  correction  of  these 
errors  by  manual  methods. 

Once  cargo  is  within  one  of  the  three  west  coast  ocean 
terminals,  it  falls  under  SURS ,  or  SURS/CARDPAC.   A  CARDPAC 
is  a  set  of  eight  IBM  cards  containing  various  shipping 
data.   These  cards  have  been  punched  by  terminal  personnel 
with  data  obtained  from  METS  and  any  additional  information 
that  may  be  needed  from  shipping  forms. 

The  CARDPAC;  the  Expected  Receipt  List  (ERL) ,  a  list 
of  the  cargo  expected  in  a  particular  shipment;  and  the 
applicable  source  document,  e.g.,  a  Government  Bill  of 
Lading  (GBL) ,  Commercial  Bill  of  Lading  (CBL) ,  Dray  Tag, 
or  Transportation  Control  and  Movement  Document  (TCMD)  are 
forwarded  to  the  terminal's  cargo  receiving  area.   As  the 
cargo  is  unloaded  a  check  is  made  against  information  con- 
tained on  the  CARDPAC  and  corrections  are  made  if  necessary. 
The  CARDPAC  is  then  attached  to  the  shipment  and  will  be 
utilized  to  monitor  its  flow  through  the  terminal  facilities 

In  summary,  by  using  the  METS  and  SURS  management  infor- 
mation systems,  WAMTMC  is  able  to  monitor  the  flow  of  cargo 
into  and  within  its  terminal  facilities.   Management  reports 
generated  by  the  system  and  statistical  analyses  based  on 
the  information  provided  by  the  system  are  only  as  good  as 
the  data  going  into  it. 
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H.   CARGO  CLASSIFICATION 

There  are  two  types  of  cargo  passing  through  the  Con- 
tainer Stuffing  Station:   (1)  release  unit  (RU) ;  and  (2) 
less  than  release  unit  (LRU) . 

Shipments  in  excess  of  10,000  pounds  and  special  cate- 
gory cargo,  such  as  classified  cargo,  are  classified  as 
release  unit  material.   Under  DoD  Regulation  4500. 32-R, 
Military  Transportation  and  Movement  Procedures  (MILSTAMP) , 
release  unit  material  requires  positive  export  traffic 
release . 

Consignors  having  break  bulk  or  containerizable  cargo 
that  fall  into  RU  classification  must  request  and  receive 
clearnace  to  move  that  material  before  it  can  be  shipped. 
WAMTMC ,  Export  Control  Division,  serves  as  the  clearance 
authority  for  RU  material  for  the  fourteen  Western  states. 

This  thesis  is  not  concerned  with  break  bulk  or  source 
stuffed  cargo.   Therefore,  this  paper  will  not  go  into  a 
detailed  description  of  the  booking  procedure  except  to 
mention  that  it  is  accomplished  through  the  coordinated  ef- 
forts of  MSCPAC,  WAMTMC,  the  shipper  and  the  commercial 
carrier. 

It  should  be  mentioned  that  if  RU  classified  cargo  is 
not  of  sufficient  quantity  to  be  source  stuffed,  it  is 
shipped  to  the  CFD  for  stuffing. 

Eighty  percent  of  the  shipments  passing  through  MOTBA 
are  of  the  LRU  classification.  LRU  material  does  not  re- 
quire positive  release  into  the  system.  Therefore,  there 
is  an  uncontrolled  flow  of  cargo  coming  into  the  Container 
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Freight  Division.   This  results  in  a  very  difficult  booking 
problem  for  the  Container  Freight  Division,  MOTBA,  WAMTMC 
and  MSCPAC.   It  is  the  forecast  of  this  uncontrolled  cargo 
flow  to  which  this  thesis  addresses  itself. 

Figure  3  represents  the  cargo  cycle  to  be  considered. 
The  procedures  at  each  step  are  as  follows  [Ref .  15] : 

1)  A  DoD  overseas  command  makes  a  request  for  goods 

by  requisition.   The  specific  requestor  is  the  user  of  these 
goods  and  the  ultimate  consignee.   The  requisition  is  to  be 
filled  by  the  shipper  (consignor)  who  may  be  a  commercial 
vendor  or  a  military  depot. 

2)  Material  that  is  classified  as  less  than  release 
unit  is  shipped  by  the  consignor  directly  to  the  Container 
Freight  Division  at  MOTBA. 

3)  Under  Military  Traffic  Management  Regulation,  con- 
tainer requirements  are  offered  to  WAMTMC  by  MOTBA. 

4)  A  request  for  booking  is  then  submitted  to  MSCPAC. 

5)  MSCPAC  then  books  space  aboard  ocean  carrier. 

6)  Once  booking  aboard  the  vessel  is  made,  the  related 
information  is  forwarded  to  WAMTMC. 

7)  The  Container  Freight  Division  at  MOTBA  receives 
the  booking  data  from  WAMTMC. 

8)  The  Container  Freight  Division  communicates  directly 
with  the  ocean  carrier  to  arrange  for  spotting  and  pick-up 
of  containers. 

9)  Once  containers  have  been  stuffed  at  the  CSS,  they 
are  delivered  directly  to  the  ocean  carrier's  container  yard, 
From  this  point  shipment  from  POE  to  POD  is  made. 
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10)   Delivery  of  material  is  made  to  the  consignee  by 
either  the  ocean  carrier  or  by  arrangements  made  by  the 
overseas  command. 

I.   CURRENT  BOOKING  PROCEDURES 

As  previously  stated,  a  request  for  booking  is  made 
twenty-one  days  in  advance  and  is  based  on  a  forecast  of 
the  volume  of  cargo  in  measurement  tons  that  is  expected  to 
be  on-hand  at  the  time  the  vans  are  to  be  stuffed.   An  of- 
fering is  made  based  on  the  following  formula: 

Booking  _  ,  Cargo-   ^    ,-  Cargo   >.  +  ,   Forecasted     ^ 
Offering    *■  on-hand  *         *•  Booked  *         *■  Cargo  Receipts  *  ' 

The  CSS  computes  the  previous  three  weeks  average  weekly 
receipts  and  uses  this  figure  as  a  basis  for  forecasting  the 
volume  of  cargo  arriving  at  the  facility  between  the  date  of 
offering  and  the  stuffing  of  containers.   Figure  3  illustrates 
the  present  method. 

Figure  4' depicts  the  situation  where  a  particular  POD 
has  100  measurement  tons  (MT's)  of  cargo  on  the  floor  on 
Julian  calendar  day  70.   To  keep  the  computations  simple, 
250  MT's  of  cargo  will  represent  the  average  weekly  receipts 
for  this  POD  for  the  past  three  weeks  (15  days,  since  only 
weekdays  are  included)  giving  average  daily  receipts  of  50 
MT's.   Container  space  has  already  been  booked  aboard  three 
vessels  during  the   21-   day  period  under  consideration. 
The  cutoff  dates  have  been  made  for  days  77,  79  and  84  for 
360,  80  and  120  MT's  respectively.   Using  the  past  three 
weeks  average  receipts  as  the  basis  of  the  forecast,  there 
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will  be  360  MT's  on-hand  on  day  77.   The  two-day  period 
where  no  cargo  accumulation  occurs  represents  a  Saturday 
and  a  Sunday.   The  five  working  days  will  accumulate  250 
MT's.   Since  space  for  360  MT's  has  been  previously  booked, 
all  cargo  on  the  floor  for  this  POD  will  be  stuffed.   Two 
days  later  there  will  be  100  MT's  on  the  floor  but  20  MT's 
will  not  be  stuffed.   On  day  84  only  120  of  the  170  MT's 
will  be  stuffed.   By  day  91,  the  day  for  which  the  forecast 
is  made,  the  forecast  says  there  will  be  300  MT's  on  the 
floor  for  this  particular  POD.   The  number  of  vans  requested 
will  depend  on  the  average  cube  utilization  for  the  POD. 
If  the  POD's  average  is  50  MT's  then  an  offering  for  six 
vans  will  be  made  on  day  70  to  be  available  for  stuffing 
on  day  91. 

J.   CANCELLATION  PROCEDURE 

Cancellation  of  space  booked  aboard  a  carrier  must  be 
made  no  later  than  "a  reasonable  length  of  time"  prior  to 
the  cutoff  date  according  to  MSC's  Container  Agreement  and 
Rate  Guide,  RG8.   If  cancellation  cannot  be  made  within  this 
time  and  if  the  carrier  cannot  utilize  this  space,  then  the 
government  may  be  charged  for  the  space  even  though  it  was 
not  used.   The  present  policy  at  WAMTMC  is  to  cancel  no 
later  than  five  days  prior  to  cutoff. 

K.   OPERATIONAL  RESTRICTIONS  AND  REQUIREMENTS 

It  should  be  noted  that  there  are  other  operational  re- 
strictions and  requirements  that  affect  stuffing  and  booking 
procedures,  although  a  thorough  understanding  of  them  is  not 
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essential  to  this  work.   Briefly  they  are:   1)  a  first-in- 
first-out  (FIFO)  cargo  stuffing  procedure  requiring  a  mini- 
mum average  FIFO  performance  of  80%  for  each  POD,  2)  shipment 
priorities  that  are  exceptions  to  the  FIFO  procedure,  3)  a 
requirement  for  the  timely  handling  of  certain  goods  such 
as  household  effects,  4)  a  requirement  for  a  monthly  cube 
utilization  average  of  75%  with  no  less  than  50%  cube  utili- 
zation to  ensure  that  container  utilization  is  more  economi- 
cal than  break  bulk  shipping,  5)  a  restriction  on  cargo 
mixing  such  as  Military  Assistance  Program  (MAP)  cargo  not 
being  permitted  to  be  mixed  with  any  other  type  of  cargo, 
6)  the  requirement  of  contracting  the  low  cost  carrier  if 
there  is  more  than  one  that  can  provide  the  required  ser- 
vice, and  7)  a  ceiling  of  4,000  measurement  tons  on  the  floor 
at  any  one  time. 
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III.   STUDY  OBJECTIVES  AND  PROCEDURES 

A.   STUDY  OBJECTIVES 

The  essence  of  this  study  is  to  observe  the  flow  of  cargo 
into  the  Container  Freight  Division  at  MOTBA  and  attempt  to 
forecast  by  statistical  and  analytical  techniques  the  expected 
cargo  inputs  at  some  future  point  in  time.   The  forecast  is 
of  the  volume  of  containerizable  cargo  passing  through  the 
CFD  destined  for  a  particular  POD.   Break  bulk  cargo  volume 
and  source  stuffed  vans  are  not  considered  in  this  analysis. 

The  necessity  for  such  a  forecast  evolves  from  the  fact 
that  the  Military  Sealift  Command  must  book  container  space 
aboard  vessels  as  much  as  three  weeks  in  advance  of  lift. 
Because  of  the  demand  for  container  shipping  in  recent  years, 
commercial  carriers  are  able  to  make  this  demand  on  their 
customers.   It  remains  to  be  seen  whether  a  shift  in  the  sup- 
ply and  demand  factors  will  alter  this  requirement  to  any 
significant  extent. 

Although  forecasting  and  forecasting  models  are  generally 
associated  with  inventory  and  production  control,  the  same 
general  principles  apply  to  the  situation  at  MOTBA.   Es- 
sentially, the  situation  is  one  of  demand  --  the  demand  by 
overseas  commands  for  material  support.   In  this  respect  the 
forecasting  problem  can  then  be  broken  down  into  components 
of  that  demand.   Examples  of  these  components  are  average 
demand,  trend  effects,  seasonal  effects,  and  noise  or  random 
effects.   By  recognizing  these  components  one  can  then  con- 
struct statistical  models  to  deal  with  them. 
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Accurate  forecasting  may  have  an  influence  on  some  of 
the  performance  variables  that  the  container  booking  system 
generates.   These  are  frequency  of  cancellations  due  to 
overbooking,  age  of  cargo  at  stuff,  cube  utilization,  and 
total  volume  of  cargo  on  the  floor. 

In  forecasting,  the  past  is  projected  into  the  future. 
A  good  forecasting  model  should  provide  information  feedback 
accurately  and  quickly  so  that  the  manager  can  make  the 
necessary  and  appropriate  adjustments.   It  should  provide  a 
means  of  reflecting  sudden  shifts  caused  by  policy  changes 
within  the  system  or  uncontrollable  exogenous  factors  which 
may  also  influence  the  system.   A  forecast  that  is  based  on 
out-dated  information  will  not  be  beneficial  to  a  manager 
who  must  book  container  space  three  weeks  in  advance. 

A  forecast  is  only  as  good  as  the  data  it  uses.   In  this 
thesis  it  is  assumed  that  the  historical  cargo  flow  data 
used  to  evaluate  the  forecasting  methods  is  accurate.   Data 
totals  obtained  from  the  computer  tapes  were  verified  against 
the  totals  that  were  manually  computed  at  MOTBA. 

B.   PROCEDURES 

As  stated  previously  the  system  variable  that  this  thesis 
addresses  is  the  volume  of  containerizable  cargo  for  a  POD 
flowing  into  the  Container  Freight  Division  at  MOTBA.   Al- 
though there  are  weight  restrictions  to  consider,  these  re- 
strictions have  not  been  significant  factors  in  the  booking 
procedure . 
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The  data  used  throughout  this  thesis  was  taken  from 
actual  operations  at  the  Container  Stuffing  Station.   The 
data  covers  an  18-month  period  dating  from  1  January  1973 
to  30  June  1974.   This  period  generated  an  enormous  volume 
of  data  which  reflected  the  magnitude  of  the  container 
stuffing  operations  at  MOTBA. 

In  order  to  handle  this  volume  of  information  the  ana- 
lysis was  done  in  four  separate  phases.   These  were: 

1)  a  distribution  audit  in  order  to  identify  the  vari- 
ables needed  in  the  analysis, 

2)  the  grouping  of  data  in  a  more  usable  format, 

3)  a  test  for  serial  correlation  in  the  data,  and 

4)  the  application  of  various  forecasting  techniques 
and  their  evaluation. 

In  performing  the  distribution  audit  some  of  the  com- 
puter programs  from  the  McCarthy  and  Carter  thesis  were 
utilized.   Only  slight  modifications  to  these  programs  were 
needed  to  manipulate  the  enormous  volume  of  data  generated 
by  the  eighteen  months  of  operations.   The  audit  identified 
the  number  of  shipments  and  the  volume  of  these  shipments. 
From  this  information  one  could  identify  major  Ports  of 
Debarkation,  thus  narrowing  the  scope  of  the  problem. 

In  order  to  arrange  the  data  into  a  more  usable  format 
a  number  of  short  computer  programs  were  written.   These 
programs  grouped  the  data  from  daily  to  weekly  totals  in 
the  aggregate  as  well  as  POD  levels.  . 

Since  this  work  dealt  with  time-series  analysis  it  was 
reasonable  to  test  the  data  for  serial  correlation  in  order 
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to  limit  the  range  of  potential  forecasting  methods.   A 
computer  program  was  utilized  to  calculate  the  statistical 
information  needed  to  evaluate  the  various  PODs . 

Several  computer  programs  were  then  used  to  apply  various 
forecasting  techniques  to  the  data.   The  simplest  techniques 
were  applied  first.   In  some  cases  these  initial  techniques 
were  expanded  to  take  into  consideration  the  components  of 
demand  previously  mentioned.   The  various  forecasts  were 
then  evaluated  based  on  the  difference  between  the  fore- 
casted and  actual  volume. 

All  the  computer  programs  used  in  the  data  analysis  were 
written  in  the  FORTRAN  IV  computer  language.   These  programs 
were  run  on  the  IBM  360/67  computer  at  the  W.  R.  Church 
Computer  Center,  Naval  Postgraduate  School,  Monterey, 
California. 
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IV.   DATA  ANALYSIS 

A.   DISTRIBUTION  AUDIT 

Raw  data  was  taken  from  a  set  of  19  computer  tapes  sup- 
plied by  WAMTMC  and  was  combined  with  existing  data  from 
previous  research  work  to  form  a  set  of  data  covering  18 
months.   The  period  covered  was  Quarters  III  and  IV  of  FY 
'73  and  all  four  quarters  of  FY  '74.   Computer  programs 
were  used  to  extract  the  pertinent  information  from  the 
tapes,  combine  it  with  existing  data  and  place  the  aggre- 
gate information  in  alphabetical  order  by  POD  and  consignee 
in  a  data  file. 

Computer  Program  Two  of  the  McCarthy  and  Carter  thesis 
was  used  to  delineate  the  data  for  further  analysis.   This 
program  determined  the  number  of  shipments,  the  total  volume, 
the  average  volume,  the  standard  deviation  of  the  volume, 
the  total  weight,  the  standard  deviation  of  the  weight,  the 
average  density  and  the  standard  deviation  of  the  density 
for  each  POD.   The  program  also  grouped  each  POD's  shipments 
into  intervals  of  volume,  weight  and  density. 

The  next  procedure  entailed  the  examination  of  the  gross 
weekly  totals  of  shipments  and  volume  of  cargo  arriving  at 
the  container  stuffing  station  during  the  18-month  period. 
The  following  information  was  obtained: 


SHIPMENTS 

Quarters  III  §  IV,  FY  '73 
FY  '74 
18-month  Totals 


111,463  shipments 
124,592  shipments 
236,055  shipments 
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VOLUME 

Quarters  III  §  IV,  FY  '73 
FY  «  74 
18-month  Totals 


8,398,931  cubic  feet 
10,070,724  cubic  feet 
18,469,655  cubic  feet 


There  was  obviously  a  tremendous  amount  of  activity  at 
MOTBA  during  the  second  half  of  FY  '73.   Over  47%  of  the 
total  shipments  and  45%  of  the  total  volume  for  the  18-month 
period  occurred  during  the  first  six  months  of  the  data  set. 

It  was  concluded  that  residual  Vietnam  conflict  effects 
were  responsible  for  generating  this  high  volume  of  activity 
during  the  first  six  months  of  the  18-month  period  and  that 
it  did  not  represent  the  present  and  the  expected  future 
flow  of  cargo  at  MOTBA.   Therefore,  it  was  decided  to  use 
only  the  data  for  FY  '74. 

The  audit  was  run  a  second  time  using  only  FY  '74  data. 
During  that  period  157  PODs  and  3,356  consignees  were 
identified  as  having  received  shipments  during  this  period. 
These  PODs  and  consignees  accounted  for  124,592  shipments 
with  a  total  volume  of  10,070,724  cubic  feet  as  previously 
stated.   This  was  an  average  weekly  total  of  2,396  shipments 
with  an  average  weekly  volume  of  193,668  cubic  feet  (or 
4,842  measurement  tons).   Total  weight  for  the  period  was 
176,213,333  pounds. 

In  reviewing  the  yearly  gross  totals  for  the  157  PODs, 
it  was  obvious  that  10%  of  these  dominated  the  activity  at 
MOTBA.  Sixteen  PODs  were  identified  as  being  the  most  ac- 
tive in  terms  of  volume  of  cargo  arriving  at  the  container 
stuffing  station  during  the  period.   The  cutoff  point  was 
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arbitrarily  chosen  at  100,000  cubic  feet  of  volume  per  year. 
As  a  result  the  least  active  of  the  16  major  PODs  averaged 
a  volume  input  of  2,041.35  cubic  feet  per  week  (or  approxi- 
mately 51  measurement  tons) .   The  number  of  shipments  for 
these  major  ports  totalled  115,732  or  92.891  of  all  ship- 
ments for  the  period.   Their  volume  total  of  9,370,693  cubic 
feet  was  93.05%  of  the  total  volume. 

Table  I,  Appendix  A,  shows  the  weekly  gross  totals  of 
shipments,  volume  and  weight  arriving  at  the  CSS  in  the  second 
half  of  FY  '73.   The  totals  for  FY  '74  appear  in  Table  II. 
Table  III,  Appendix  A,  shows  the  data  summary  of  the  dis- 
tribution audit.   Tabular  data  for  the  16  major  PODs  is 
presented  in  Table  IV.   Table  V  lists  the  activities  of  the 
major  PODs  by  shipments,  volume  and  weight. 

The  computer  output  of  the  intervals  of  weight  into 
which  each  of  the  shipments  fell  provided  a  look  at  the  por- 
tion of  the  release  unit  shipments  arriving  into  the  system. 
If  the  number  of  release  units  were  significantly  large,  it 
may  be  possible  to  forecast  with  greater  accuracy  the  volume 
of  cargo  expected  to  arrive  during  a  specified  period. 

There  were  1,900  shipments  of  ten  thousand  pounds  or 
more  during  FY  '74  using  only  data  generated  by  the  16  major 
PODs.   These  shipments  should  have  been  classified  as  re- 
lease unit  shipments  and  it  is  assumed  that  they  were.   They 
represented  1.64%  of  the  shipment  totals  for  the  major  ports. 

It  was  concluded  that  the  amount  of  release  unit  ship- 
ments arriving  at  the  container  stuffing  station  was  insignifi- 
cant in  view  of  these  figures.   Furthermore,  operations  at 
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MOTBA  indicate  that  the  expected  arrival  date  of  RU  cargo 
was  an  unreliable  item.   Therefore,  RU  cargo  was  considered 
part  of  the  uncontrolled  cargo  flow  and  grouped  with  LRU 
data  for  the  purpose  of  this  work.   Table  VI,  Appendix  A, 
shows  the  number  of  RU  shipments  classified  because  of  weight 
for  each  of  the  major  PODs . 

It  should  be  noted  that  no  attempt  was  made  to  determine 
the  number  of  shipments  classified  as  RU  for  reasons  other 
than  weight  (i.e.,  classified  material).   It  was  assumed 
that  the  volume  of  this  cargo  was  of  insignificant  proportion 
when  compared  to  the  total  yearly  volume. 

In  summary,  the  distribution  audit  provided  a  means  of 
taking  raw  data  and  arranging  it  into  a  meaningful  and  use- 
ful format.   It  permitted  the  identification  of  the  most 
active  ports  of  debarkation  which  then  became  the  focal 
point  for  further  analysis.   The  grouping  of  the  data  into 
weekly  intervals  by  POD  indicated  the  tremendous  volume 
fluctuations  that  occurred  throughout  the  period  and  the  ap- 
parent complexities  that  arise  when  attempting  to  forecast 
future  volume  flow. 

Up  to  this  point,  only  the  simple  computer  manipulation 
of  data  had  been  accomplished  through  the  distribution  audit. 
With  this  phase  completed,  there  were  still  many  questions 
to  be  answered.   What  type  of  data  had  been  generated?   Was 
there  any  correlation  in  the  volume  flow  from  week  to  week? 

The  above  questions  are  answered  in  the  following  section 
which  begins  the  statistical  analysis  of  the  data.   The  ans- 
wers then  provide  statistical  justification  for  applying 
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certain  forecasting  techniques  which  are  covered  in  section 
C  of  this  chapter. 

B.   SERIAL  CORRELATION  TEST 

The  flow  of  uncontrolled  cargo  volume  into  the  Container 
Freight  Division  at  MOTBA  is  time-series  data.   It  is  not 
uncommon  for  future  observations  in  real  time-series  data  to 
be  dependent  upon  current  observations.   This  dependence  may 
be  measured  by  analysis  of  the  consecutive  disturbances  around 
a  simple  least-square-regression  line.   If  these  disturbances 
are  found  to  be  correlated  from  one  observation  to  the  next, 
and  thus  not  random,  serial  correlation  is  said  to  be  pre- 
sent.  If  the  presence  of  serial  correlation  can  be  properly 
identified  it   can  result  in  a  more  efficient  estimation 
process  in  comparison  to  simple  least-squares  estimators. 
Such  a  process  would  give  more  weight  to  recent  data  for 
creating  predictions  of  the  values  of  future  observations. 

Johnston  [Ref.  9]  lists  three  main  consequences  that 
result  from  applying  straightforward  least-squares  formulas 
directly  to  observations  that  contain  serially  correlated 
disturbances.   They  are: 

1)  The  sampling  variances  of  the  intercept  and  the  slope 
coefficient  may  be  unnecessarily  large  in  comparison  to  those 
obtained  by  other  methods. 

2)  The  usual  least-squares  formulas  for  the  sampling 
variances  of  the  regression  coefficients  are  no  longer  valid 
since  a  serious  underestimate  of  the  variance  is  likely  to 
occur. 
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3)   The  sampling  variances  of  the  prediction  will  be 
unnecessarily  large. 

1.   The  Durbin-Watson  d-Statistic 

When  making  forecasts  using  time-series  data  it 
would  be  a  serious  error  to  assume  serial  independence  of 
the  disturbance  term.   In  order  to  avoid  such  an  error,  the 
Durbin-Watson  d-Statistic  [Ref.  6]  provided  a  suitable  test 
for  the  presence  of  serially  correlated  disturbances.   The 
Durbin-Watson  d-Statistic  is  defined  by 


_  ti2  <zt  -  zt-i)2 


Si  fv2 


where  Z. (t=l , 2 , . . . ,n)  denotes  the  residuals  from  a  fitted- 
leas t-squares  regression.   To  determine  positive  serial 
correlation  the  observed  value  of  d  is  compared  against 
lower  and  upper  bounds  of  dT  and  d..  that  are  tabulated  for 
various  values  of  n  (the  number  of  observations)  and  k  (the 
number  of  explanatory  variables) .   The  d-statistic  for  nega- 
tive serial  correlation  is  calculated  by  taking  4  -  d. 
This  figure  is  also  compared  with  d.  and  d„.   Since  the  sign 
of  the  serial  correlation  was  unknown,  two-sided  tests  were 
conducted  by  combining  single-tail  tests.   The  test  has  four 
possible  outcomes:   (1)  positive  serial  correlation  of  the 
data;  (2)  negative  serial  correlation  of  the  data;  (3)  in- 
conclusive evidence  of  serial  correlation  with  further 
observations  being  ideally  required;  or  (4)  no  presence  of 
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serial  correlation  (which  implies  that  the  disturbances  are 
random) . 

The  calculations  were  made  using  three  different 
intervals.   The  first  interval  used  weekly  data  resulting 
in  52  observations.   The  second  interval  used  biweekly 
totals  beginning  with  the  observation  for  the  first  two 
weeks  of  FY  '74.   The  total  number  of  observations  in  this 
set  was  26.   Biweekly  totals  beginning  with  the  observation 
of  the  second  and  third  weeks  of  FY  '74  (i.e.,  ignoring  the 
first  week's  volume  total)  was  used  for  the  third  interval. 
Thus  there  were  only  25  observations  resulting  from  this 
grouping  scheme. 

The  following  example  may  help  to  clarify  the  above 
procedure.   Given  that  the  weekly  volume  input  for  POD,  BA3, 
for  the  first  seven  weeks  of  FY  '74  is  10,  20,  30,  40,  50, 
60  and  70  cubic  feet  respectively,  the  first  seven  observa- 
tions for  the  weekly  interval  data  set  correspond  according- 
ly.  The  first  three  observations  for  the  first  biweekly 
data  set  will  be  30,  70  and  110  cubic  feet.   However,  the 
first  three  observations  for  the  second  biweekly  data  set 
will  be  50,  90  and  130  cubic  feet. 

The  d-statistic  test  was  implemented  through  Com- 
puter Program  One.   In  addition  to  displaying  the  results 
of  the  calculations,  the  program  displays  the  meaning  of 
the  d-statistic  (i.e.,  positive,  negative,  inconclusive  or 
none).   Computer  Output  One  at  the  end  of  this  thesis  is  a 
sample  of  the  printout  for  Computer  Program  One. 
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2 .   Results  of  Serial  Correlation  Test 

While  the  Durbin-Watson  d-Statistic  was  calculated 
for  all  157  PODs ,  the  d- statistic  ■  £01  the  16  most  active 
PODs  was  the  focal  point  for  the  analysis. 

The  results  of  the  test  for  the  weekly  interval  data 
set  having  52  observations  for  each  POD  indicated  positive 
serial  correlation  in  the  weekly  volume  inputs  for  nine  PODs. 
The  test  was  inclusive  for  four  PODs  and  indicated  no  serial 
correlation  for  three  others.   There  were  no  PODs  with  nega- 
tive serial  correlation.   In  fact,  there  was  no  indication 
of  negative  serial  correlation  under  any  of  the  circumstances 
tested  for  all  157  PODs.   The  following  table  shows  the  test 
results  using  only  the  16  major  PODs.   The  test  was  conducted 
at  a  level  of  significance  of  0.05. 

Week  per   Week  Results  No 

Interval    Started   Positive   Negative   Incon   Correlation 

1  19         0        4        3 

2  15         0        5        6 
2  2         7         0        5        4 

Table  VII  gives  the  d-statistic  results  for  each  POD 
for  the  three  intervals  tested. 

Given  the  results  of  the  d-statistic  tests  it  was 
then  decided  to  test  for  the  aggregate  presence  of  serial 
correlation  in  the  data  for  the  16  major  PODs.   The  null 
hypothesis  for  the  test  was 

HQ:  Pl  =  p2  =  ...  =  P16  =  0 
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where  p-,...,p1fi  symbolically  represents  the  serial  correla- 
tion coefficients  for  the  16  major  PODs .   The  alternate 
hypothesis  was 

H-:   not  all  p's  =  0. 

This  hypothesis  is,  in  essence,  testing  the  plausi- 
bility that  the  individual  cases  indicating  positive  serial 
correlation  were  Type  I  errors,  i.e.,  a  rejection  of  the 
null  hypothesis,  p.  =0,  when  it  was  true.   This  was  formu- 
lated as  a  Bernoulli  process  with  the  probability  of  success 
(the  indication  of  positive  serial  correlation)  equal  to  the 
level  of  significance  for  the  d-statistic  test  (0.05). 

Since  a  Bernoulli  process  was  involved,  the  binomial 
distribution  function  was  applied.   Letting  the  probability 
of  success,  p,  equal  0.05;  the  number  of  trials,  n,  equal 
to  16;  and  the  specific  number  of  successes,  r,  equal  to  5, 
7  and  9  respectively,  the  decision  rule  was 

Accept  H   if  P(R>r|n,p)  >  0.05 

Reject  H   if  P(R>r|n,p)  <  0.05 

where  R  equals  the  unknown  number  of  successes. 

Using  cumulative  probability  distribution  tables  the 
following  results  were  obtained: 

P(R>5| 16,0.05)  =  0.0009 
P(R>7|16,0.05)  =  0.0001 
P(R>9| 16,0.05)  <  0.0001. 

Based  on  these  results  the  null  hypothesis  was  re- 
jected and  it  was  assumed  that  significant  serial  correlation 
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surely  existed  somewhere  in  the  data.   Therefore,  there  was 
justification  for  using  methods  other  than  ordinary-least- 
squares  regression  to  forecast  the  volume  of  cargo  arriving 
at  the  Container  Freight  Division  at  MOTBA. 

C.   FORECASTING  METHODS 

The  results  of  the  Durbin-Watson  d-Statistic  indicated 
the  presence  of  serial  correlation  and  thus  the  serial  de- 
pendence of  the  observations.   It  showed  that  the  distur- 
bances around  an  ordinary-least-squares  regression  line 
would  not  fall  equally  on  each  side  of  the  line.   As  the 
disturbances  tended  to  fall  more  on  one  side  of  the  line 
than  the  other,  this  implied  that  the  most  recent  distur- 
bance was  the  best  predictor  of  the  subsequent  distrubance. 
Also,  by  implication,  the  best  predictor  of  the  next  obser- 
vation is  the  most  recent  observation  or  observations.   This 
supplied  justification  for  the  application  of  forecasting 
methods  other  than  ordinary- least-squares  to  place  more 
emphasis  on  recent  data  rather  than  the  minimization  of  the 
squared  residual  term.   While  there  were  a  number  of  methods 
from  which  to  choose,  the  methods  used  in  this  work  were 
chosen  because  of  their  relevency  to  the  problem  and  for 
their  ease  of  calculation. 

This  section  looks  at:   the  moving  average  method  that 

is  presently  being  used  at  MOTBA;  the  exponentially  weighted 

average  method;  the  exponentially  weighted  average  method 

with  adjustment  for  trend;  and  an  adaptive  response  rate 

technique  applied  to  the  exponentially  weighted  average 

method. 
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These  methods  were  evaluated  using  mean  absolute  devia- 
tion (MAD)  and  mean  deviation  squared  (MDS)  as  the  criteria 
for  goodness.   By  computing  the  mean  of  the  deviations  it 
was  possible  to  compare  the  various  forecasting  methods  when 
they  resulted  in  differing  numbers  of  forecasts  generated 
during  the  period. 

In  addition,  MAD  and  MDS  were  computed  for  one-,  two- 
and  three-week  time  horizons.   For  the  one-week  time  horizon 
the  forecast  was  subtracted  from  the  actual  input  volume  on 
a  week  by  week  basis.   For  the  two- (three- )week  time  horizon 
the  latest  forecast  was  multiplied  by  two  (three)  and  sub- 
tracted from  the  actual  input  volume  for  the  two  (three) 
weeks  covered  by  the  forecast.   This  was  done  on  a  week  by 
week  basis. 

The  following  example  may  help  to  clarify  time  horizon 
calculations.   Say  that  the  forecast  for  week  20  is  1,000 
cubic  feet.   If  a  three-week  time  horizon  is  being  used,  the 
forecast  will  be  3,000  cubic  feet  for  weeks  20,  21  and  22. 
If  the  actual  input  volumes  turn  out  to  be  900;  1,100  and 
1,200  cubic  feet  respectively,  the  absolute  deviation  for 
the  forecast  is  then  200  cubic  feet  (or  900  +  1,000  +  1,200 
-  3,000).   Week  21  will  result  in  a  new  forecast  for  weeks 
21,  22  and  23.   The  forecast  for  week  21  is  multiplied  by 
three  and  subtracted  from  the  actual  input  volumes  for  weeks 
21,  22  and  23  to  determine  the  absolute  deviation  for  that 
particular  forecast.   Thus,  for  each  forecast  there  is  an 
absolute  deviation  that  is  computed. 
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What,  then,  was  the  purpose  for  using  time  horizons? 
One  reason  was  the  requirement  for  booking  container  space 
three  weeks  in  advance  and  thus  projecting  a  one-week  fore- 
cast over  a  period  greater  than  one  week.   One-,  two-  and 
three-week  horizons  were  used  to  make  comparisons  of  any 
effects  that  may  be  caused  by  changes  in  the  booking  re- 
quirements.  It  was  also  anticipated  that  an  increase  in 
the  time  horizon  may  have  some  effect  on  the  results  of  the 
various  forecasting  methods  due  to  a  dampening  effect  on 
data  with  large  fluctuations.   This  dampening  effect  is  ex- 
plained later. 

It  should  also  be  noted  that  comparisons  among  the  var- 
ious forecasting  methods  were  only  made  within  a  given  time 
horizon. 

Tables  VIII  and  IX,  Appendix  A,  show  the  MADs  and  MDSs 
generated  by  all  of  the  forecasting  methods  used  in  this 
analysis . 

1 .   Moving  Average  Method 

As  previously  explained,  the  CSS  at  MOTBA  computes 
the  previous  three  weeks  average  receipts  in  order  to  fore- 
cast the  input  volume  that  will  be  on  hand  in  21  days.   The 
offerings  for  bookings  that  are  made  depend  on  three  factors 
1)  the  volume  of  cargo  on  hand;  2)  the  volume  of  cargo 
booked;  and  3)  the  forecasted  volume. 

The  three-week  moving  average  is  nothing  more  than 
a  linear  average.   This  average,  multiplied  by  three  is  then 
used  to  forecast  the  volume  receipts  for  the  next  three 
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weeks.   In  this  method  equal  weight  is  given  to  each  of  the 
past  three  weekly  volume  totals. 

In  order  to  evaluate  the  moving  average  method  the 
MADs  and  MDSs  were  calculated  for  one-week  to  ten-week 
moving  averages  for  the  16  PODs .   In  addition,  the  aggre- 
gate MAD  and  MDS  for  each  POD  were  calculated  for  one-  to 
ten-week  moving  averages.   In  other  words,  a  specific  moving 
average  (from  one  to  ten  weeks)  was  applied  to  the  volume 
inputs  for  each  POD.   The  MAD  and  MDS  for  each  POD  were 
calculated  and  then  summed  to  get  the  aggregate.   The  same 
calculations  were  made  for  two-  and  three-week  time  horizons 

Computer  Program  Two  was  utilized  to  obtain  the 
calculations  for  the  analysis  of  the  moving  average  method. 
Computer  Output  Two  shows  a  portion  of  the  output  that  was 
generated  by  the  program. 

a.   Data  Analysis  of  the  Moving  Average  Method 

In  13  of  16  cases  a  moving  average  other  than 
the  three-week  moving  average  provided  a  better  forecast 
when  MAD  was  the  criterion  of  goodness.   When  using  MDS  as 
the  criterion,  the  result  was  15  of  16  cases.   The  time  in- 
tervals that  provided  the  best  forecast  will  be  called  the 
"optimal  moving  average." 

These  evaluations  were  made  using  a  one-week 
time  horizon.   Similar  results  were  obtained  using  two-  and 
three-week  time  horizons  except  that  the  optimal  time  inter- 
val for  the  moving  average  was  not  the  same  in  each  case. 

It  was  noted  that  in  some  cases  the  selection  of 
the  goodness  criterion  determined  the  optimal  interval  for 
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the  moving  average.   This  was  the  case  for  six  PODs  using 
a  one-week  time  horizon  and  was  true  to  a  greater  or  lesser 
extent  for  the  two-  and  three-week  time  horizons. 

The  decrease  in  the  MAD  when  a  moving  average 
other  than  a  three-week  moving  average  was  optimal,  varied 
from  4.5%  to  18.5%.   Similar  reductions  in  MAD  and  MDS  oc- 
curred when  the  other  time  horizons  were  used. 

The  optimal  time  interval  for  a  moving  average 
that  applied  to  all  PODs  was  seven,  ten  and  ten  weeks  for 
one-,  two-  and  three-week  time  horizons,  respectively,  with 
MAD  as  the  criterion.   The  same  time  intervals  held  true 
with  MDS  as  the  criterion. 

It  was  noted  that  the  number  of  forecasts  de- 
creased as  the  time  interval  for  the  moving  averages  in- 
creased.  For  example,  in  a  three-week  moving  average,  the 
first  three  weeks  are  averaged  to  make  a  forecast  for  week 
four.   Thus,  there  were  only  49  forecasts  made  for  the 
period.   For  a  ten-week  moving  average  there  were  only  42 
forecasts.   The  use  of  two-  and  three-week  time  horizons 
resulted  in  even  fewer  forecasts.   Comparison  among  the 
various  moving  averages  for  a  specific  time  horizon  was 
possible  since  mean  absolute  deviation  and  mean  deviation 
squared  were  used  as  the  criteria  for  goodness. 

As  the  time  interval  for  the  moving  average  in- 
creased, the  greater  was  the  chance  for  large  fluctuations 
in  the  data  to  be  dampened  by  the  initial  forecast.   For 
example,  the  volume  inputs  for  the  first  six  weeks  of  FY 
'74  for  TA1  were  2,611;  4,307;  1,079;  3,360;  6,765  and 

56 


2,429  cubic  feet,  respectively.   These  fluctuations  were 
dampened  when  moving  averages  greater  than  three  weeks  were 
applied  to  the  data  (with  a  five-week  average  being  opti- 
mal) .   This  dampening  effect  stressed  the  fact  that  a  three- 
week  moving  average  placed  more  emphasis  on  recent  data. 
Longer  time  intervals  placed  less  emphasis  on  recent  data 
and  considered  older  data  as  well. 

It  was  also  noted  that  the  MAD  and  MDS  decreased 
as  the  time  horizon  increased.   The  reason  was  that  the 
two-  and  three-week  time  horizons  also  tended  to  cancel 
the  volume  fluctuations.   A  look,  again,  at  Figure  1  shows 
the  very  large  volume  fluctuations  that  occur  in  the  POD 
data  and  how  it  is  possible  for  the  dampening  effect  to  take 
place. 

2 .   Exponentially  Weighted  Average  Method 

In  lieu  of  the  results  of  the  serial  correlation 
test,  a  viable  alternative  to  the  moving  average  method  is 
the  exponentially  weighted  average  (EWA) .   In  its  simplest 
form  an  EWA  is  based  on  a  period  by  period  adjustment  of  the 
most  recent  forecasted  average  by  adding  or  subtracting  a 
fraction  of  the  difference  between  the  actual  volume  in  the 
current  period  and  the  latest  forecasted  average  volume. 
It  differs  from  the  moving  average  in  that  the  weight  given 
the  most  recent  observation  can  be  varied  while  also  giving 
consideration  to  all  past  data.   EWA  is  represented  by 


t 


=  a  V   +  (1  -  a)  F 
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where  F   is  the  forecast  average,  V\  is  the  actual  input 
volume  and  a  is  the  exponential  smoothing  constant  which 
must  be  between  zero  and  one.   The  forecast  for  the  next 
period  F   -  can  be  taken  directly  from  the  computed  value 
of  F  .   This  can  be  justified  because  trends  and  seasonal 
adjustments  are  not  accounted  for  in  the  model  [Ref .  7] . 

The  exponentially  weighted  average  method  is  also 
known  as  the  exponential  smoothing  method.   Both  names  are 
used  in  this  thesis. 

The  calculations  for  the  EWA  method  were  implemented 
by  Computer  Program  Three.   The  MAD  and  MDS  for  each  POD 
were  calculated  for  a-levels  between  -0.10  and  +1.00  in  in- 
crements of  0.01.   In  addition  the  MAD  and  MDS  were  computed 
for  a-levels  that  were  applied  to  all  16  PODs.   In  other 
words  a  specific  a-level  was  applied  to  data  for  each  POD. 
The  MAD  and  MDS  for  each  POD  were  calculated  and  then  summed 
to  obtain  an  aggregate  total.   The  a-levels  which  generated 
the  smallest  MAD  and  MDS  are  referred  to  as  the  "optimal 
a-levels  for  aggregate  data."   It  should  be  noted  that  the 
optimal  a-levels  for  aggregate  data  may  differ  depending  on 
the  criterion  of  goodness. 

A  canned  computer  subroutine,  SUBROUTINE  OSPLOT, 
available  at  the  computer  center,  was  used  to  plot  the 
a-levels  versus  the  MAD  for  each  POD.   The  plot  displayed 
how  the  MAD  reacted  to  various  levels  of  a.   Thus  it  was 
possible  to  visually  determine  the  optimal  a-level  for 
aggregate  data  with  MAD  as  the  criterion  of  goodness. 
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Computer  Output  Three -Two  shows  the  MADs  and  MDSs 
calculated  for  various  a-levels.   Computer  Output  Three-One 
is  a  sample  of  the  OSPLOT  for  four  of  the  PODs . 
a.   Data  Analysis  of  EWA  Method 

With  MAD  as  the  criterion  for  goodness,  the  EWA 
method  resulted  in  a  decrease  in  MAD  in  14  of  16  cases  when 
compared  to  the  three-week  moving  average.   The  decrease 
ranged  from  4.9%  for  RJ3  to  12.41  for  XJ1.   However,  when 
EWA  was  compared  to  the  optimal  moving  averages  for  the 
individual  PODs  it  produced  a  smaller  MAD  in  only  4  of  16 
cases.   The  a-levels  varied  from  0.04  to  0.60  for  the  16 
PODs. 

The  optimal  a-level  for  aggregate  MAD  data  was 
0.36.   The  aggregate  MAD  for  the  EWA  method  was  59,387 
cubic  feet  as  compared  to  60,172  cubic  feet  for  the  three- 
week  moving  average  and  59,209  cubic  feet  for  the  optimal 
moving  average. 

The  preceding  analysis  was  made  from  calcula- 
tions using  a  one-week  time  horizon.   Similar  results  were 
obtained  using  two-  and  three-week  horizons.   It  is  interest 
ing  to  note  that  there  were  changes  in  the  a-levels  for  most 
PODs  depending  on  the  time  horizon.   One  POD  (UM4)  produced 
negative  a-levels  (-0.02  and  -0.01)  for  two-  and  three- 
week  time  horizons,  respectively.   The  optimal  a-levels  for 
aggregate  data  were  0.35  for  the  two-week  horizon  and  0.33 
for  the  three-week  horizon.   As  was  the  case  with  the  moving 
average  method,  the  aggregate  MADs  decreased  as  the  time 
horizons  increased  due  to  the  dampening  of  the  fluctuations. 
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With  MDS  as  the  criterion  for  goodness  the  EWA 
method  resulted  in  a  decrease  in  MDS  in  12  of  16  cases  when 
compared  to  the  three-week  moving  averages.   When  compared 
to  the  optimal  moving  average,  there  was  a  decrease  in  only 
5  of  16  cases.   In  each  case  the  optimal  a-level  for  MDS 
was  different  from  the  optimal  level  for  MAD;  although  the 
optimal  a-level  for  aggregate  data  was  approximately  equal 
for  the  MDS  (0.37)  and  the  MAD  (0.36).   The  aggregate  MDS 
for  the  EWA  method  was  less  than  the  MDS  for  the  three-week 
moving  average  but  greater  than  the  MDS  for  the  optimal 
moving  average. 

A  one-week  time  horizon  was  used  in  the  above 
analysis  for  MDS.   As  was  the  case  with  the  MAD  analysis 
of  individual  PODs ,  two-  and  three-week  time  horizons 
generated  similar  results.   However,  for  a  two-week  time 
horizon  of  aggregate  data,  EWA  generated  a  slightly  smaller 
MDS  than  that  generated  by  the  optimal  moving  average. 
This  was  not  the  case  for  the  three-week  time  horizon, 
however. 

The  way  in  which  the  EWA  method  was  applied  re- 
sulted in  no  forecast  deviations  being  computed  until  the 
third  week.   This  resulted  in  some  smoothing  of  the  data  but 
not  as  much  smoothing  as  occurred  with  moving  averages  with 
intervals  greater  than  three  weeks.   When  applying  the  EWA 
method  to  data  with  large  fluctuations  at  the  beginning  of 
the  year,  there  was  a  lag  in  the  forecast  until  these  fluc- 
tuations were  dampened  over  time.   This  was  the  apparent 
reason  for  the  seven-week  moving  average  being  a  better 

60 


forecaster  than  the  EWA  for  a  one-week  time  horizon  and  a 
ten-week  moving  average  being  better  for  a  two-  and  three- 
week  time  horizon. 

In  summary,  the  EWA  method  provided  a  means  of 
emphasizing  past  data  as  well  as  placing  some  significance 
on  the  most  recent  data  by  the  use  of  the  smoothing  con- 
stant.  The  EWA  method  provided  a  better  forecast  than  the 
three-week  moving  average;  however,  the  decrease  in  the  MAD 
for  aggregate  data  was  only  2.2%.   For  individual  PODs  this 
decrease  was  as  high  as  18%,  but  for  other  PODs  there  was 
an  increase  in  the  MAD  when  EWA  was  used. 

3.   Exponentially  Weighted  Average  Method  with  Trend 
Adjustment 

Neither  the  moving  average  or  the  EWA  took  into 
consideration  any  adjustment  for  trends  that  may  exist  in 
the  data.   Also,  the  slopes  of  the  least-squares -regression 
line  for  some  PODs  used  in  the  d-statistic  calculations  gave 
an  indication  of  the  possibility  of  trend  effects  in  the 
data.   Therefore,  the  EWA  was  modified  to  take  into  account 
the  adjustment  for  trend. 

The  difference  between  the  volume  forecast  averages 
from  week  to  week  is  the  apparent  trend.   Exponential  smooth 
ing  can  be  applied  to  this  trend  just  as  it  was  applied  to 
the  forecast  average.   This  trend  adjustment  is  represented 
by 
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By  combining  this  trend  adjustment  with  the  new  forecast 
average,  the  forecast  for  the  next  period  is  then 
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where  F   -.  is  the  forecast  for  the  exponentially  weighted 
average  with  trend  adjustment  (EWA-TA)  [Ref.  7], 

Computer  Program  Four  was  used  to  make  the  statisti- 
cal calculations  for  the  EWA-TA.   Except  for  the  trend  ad- 
justment, the  computational  procedures  were  the  same  as  those 
for  the  EWA  method.   Once  again,  SUBROUTINE  OSPLOT  was  used 
to  plot  a-levels  versus  MAD  for  each  POD. 

a.   Data  Analysis  of  EWA-TA 

Using  a  one-week  time  horizon  and  MAD  as  the 
criterion,  a  decrease  in  MAD  resulted  in  only  3  of  16  cases 
when  comparing  EWA-TA  to  EWA.   Only  in  the  case  of  RG1  was 
there  any  significant  change  (an  8%  decrease).   However, 
this  POD  was  an  exception  to  the  rule  of  fluctuating  volume 
inputs.   Its  weekly  volume  inputs  toward  the  last  half  of 
the  year  had  remained  within  a  relatively  small  range.   In 
only  one  case  did  EWA-TA  provide  a  better  forecast  than  any 
of  the  methods  previously  discussed.   However,  the  improve- 
ment was  minimal.   In  half  the  cases  the  difference  in  MADs 
between  EWA-TA  and  EWA  was  less  than  50  cubic  feet.   Also, 
the  a-level  dropped  significantly  in  all  16  cases.   The  op- 
timal a-level  for  aggregate  data  was  0.16;  however,  the  MAD 
was  greater  by  1,904  cubic  feet. 
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With  MDS  as  the  criterion  for  goodness  the  re- 
sults were  similar.   Improvement  on  the  EWA  method  occurred 
in  only  two  cases.   In  most  cases  the  a-levels  were  differ- 
ent from  those  with  MAD;  however,  the  optimal  a-level  for 
aggregate  data  changed  by  only  0.01  to  0.17.   The  aggregate 
MDS  was  the  largest  of  all  the  methods  applied  thus  far. 
This  was  true  for  all  time  horizons. 

Very  few  changes  occurred  when  two-  and  three- 
week  time  horizons  were  applied.   Once  again,  UM4  produced 
negative  a-levels.   The  optimum  a-levels  for  aggregate  data 
did  not  change  significantly.   The  aggregate  MAD  and  MDS 
decreased  in  the  same  proportion  as  the  aggregate  totals  for 
the  other  methods. 

In  summary,  the  exponentially  weighted  average 

method  with  trend  adjustments  did  not  provide  an  improved 

forecast  over  the  moving  average  method  or  the  exponentially 

weighted  method  without  trend  adjustments.   Apparently  there 

was  no  definite  trend  in  the  volume  inputs  for  the  PODs .   It 

is  doubtful  that  the  method  can  be  effective  when  applied 

to  data  covering  only  a  one  year  period. 

4 .   Exponentially  Weighted  Average  with  an  Adaptive 
Response  Rate 

Because  of  the  large  volume  fluctuations  in  the  data, 

the  application  of  an  adaptive  response  rate  (ARR)  provided 

a  plausible  alternative  to  a  simple  exponentially  weighted 

average.   The  ARR  provides  an  automatic  means  of  detecting 

sudden  volume  changes  through  the  use  of  a  tracking  signal. 

Trigg  and  Leach  [Ref.  18]  define  the  tracking  signal  as 
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t~„„i,-    c  •    i       Smoothed  Error 

Tracking  Signal  =  * .u    A    .« = — - — E 

to   to      Smoothed  Absolute  Error 

where     the  error  is  the  difference  between  the  forecasted 
volume  and  the  actual  volume  inputs  in  each  period.   The 
error  was  smoothed  by  choosing  a  value  for  gamma,  y,    between 
-1  and  +1  such  that 

New  Smoothed  Error  =  (1-y)  old  smoothed  error  +  y   error, 

New  Smoothed  Absolute    f~      ^      -,  ,      .,  ,   ,   , 

r.  =  (1-y)  old  smoothed  absolute  error 

Error  v   ■  ' 

+  y  absolute  error. 

The  tracking  signal  was  then  applied  to  the  EWA  methods  with- 
out trend  adjustment  by  letting  the  smoothing  constant  equal 
the  absolute  value  of  the  tracking  signal.   Thus,  it  ap- 
peared that  ARR  would  provide  a  rapid  response  to  large 
volume  changes  whereas  the  constant  model  would  lag  behind 
these  sudden  changes. 

The  adaptive  response  rate  was  also  applied  to  the 
exponentially  weighted  average  with  trend  adjustment.   It 
resulted  in  very  large  deviations  because  of  the  conflict 
between  the  trend  adjustment  and  the  response  rate.   They 
were  apparently  working  against  each  other;  and  therefore, 
this  method  will  be  ignored  in  this  paper. 

MADs  and  MDSs  for  the  16  major  PODs  were  calculated 
for  ylevels  between  -1.00  and  +1.00  in  increments  of  0.01. 
Optimal  ylovels  for  aggregate  data  were  also  computed. 
These  calculations  were  implemented  by  Computer  Program 
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Five.   Plots  of  ylevels  versus  MAD  were  made  through  the 
use  of  the  special  subroutine 

a.   Data  Analysis  of  EWA  with  ARR 

With  a  one-week  time  horizon  and  MAD  as  the 
criterion, the  ARR  resulted  in  the  smallest  MAD  in  6  of  16 
cases  when  compared  to  all  other  methods  previously  dis- 
cussed.  The  reduction  in  MAD  was  less  than  5%  in  five  of 
these  cases.   The  optimal  y-level  for  aggregate  data  was 
-0.05.   The  aggregate  MAD  at  this  level  was  smaller  than 
the  MADs  for  the  three-week  moving  average  and  the  EWA-TA; 
but  it  was  larger  than  the  MADs  for  the  optimal  moving 
average  and  the  EWA  method. 

A  two-week  time  horizon  produced  similar  re- 
sults on  the  individual  POD  as  well  as  aggregate  level. 
Individual  POD  results  were  similar  with  a  three-week  time 
horizon.   For  aggregate  data  only  the  optimal  moving  aver- 
age resulted  in  a  smaller  MAD.   However,  the  reduction  was 
less  than  II,  and  for  all  practical  purposes  it  was  equal. 

With  a  one-week  time  horizon  and  MDS  as  the 
criterion,  the  ARR  produced  a  lower  MDS  in  5  of  16  PODs 
when  compared  to  all  other  methods.   Two-  and  three-week 
time  horizons  had  similar  results  on  the  POD  level.   For 
aggregate  data,  the  ARR  was  unable  to  generate  an  MDS  smaller 
than  any  of  the  other  methods  for  the  three  time  horizons. 

In  summary,  the  ARR  provided  a  more  complicated 
approach  to  exponential  smoothing  than  any  of  the  other 
methods.   On  the  individual  POD  level  it  provided  a  better 
forecast  in  some  cases;  yet,  the  improvement  was  minimal. 
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The  time  horizon  and  the  criterion  for  goodness  were  signif- 
icant factors  in  determining  whether  the  ARR  generated  the 
better  or  best  forecast. 

While  the  purpose  of  the  ARR  was  to  respond  more 
rapidly  to  large  fluctuations  in  the  data,  it  would  be  more 
appropriate  to  state  that  it  was  designed  to  respond  to 
large  fluctuations  which  resulted  in  a  shift  in  the  data. 
While  there  were  large  fluctuations  in  the  data  there  were 
no  apparent  shifts  in  the  data.   Plots  of  the  weekly  volume 
inputs  indicated  that  there  were  large  fluctuations  in  the 
data,  but  that  these  fluctuations  did  not  result  in  any 
significant  shifts  in  the  data. 

5.   Summary  of  Forecasting  Methods 

Four  forecasting  techniques  were  applied  to  the 
data:   1)  the  moving  average  method;  2)  the  exponentially 
weighted  average  method;  3)  the  exponentially  weighted 
average  method  with  trend  adjustment;  and  4)  the  exponen- 
tially weighted  average  method  with  an  adaptive  response 
rate.   Past  data  was  utilized  to  forecast  future  inputs  in 
each  of  the  methods.   Each  method  provided  a  means  of 
varying  the  amount  of  emphasis  placed  on  the  most  recent 
data,  forecast  average,  or  error. 

In  order  to  evaluate  the  forecasting  methods,  the 
mean  absolute  deviation  and  mean  deviation  squared  were 
calculated  and  used  as  the  criterion  of  goodness.   The  mean 
of  the  deviations  was  used  since  each  method  resulted  in 
a  different  number  of  forecasts  being  made  for  each  period. 
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The  MADs  and  MDSs  that  were  generated  by  the  various 
forecasting  methods  indicated  the  uniqueness  of  the  volume 
inputs  for  each  POD.   The  method  providing  the  best  fore- 
cast for  a  POD  depended  upon  the  criterion  of  goodness  and 
the  time  horizon. 

Tables  VIII,  IX,  X  and  XI,  Appendix  A,  are  a  pre- 
sentation of  the  data  that  has  been  previously  discussed. 
Table  VIII  gives  the  MAD  of  the  various  methods  for  the  16 
major  PODs  in  one-,  two-  and  three-week  time  horizons. 
Table  IX  gives  the  MDS  in  the  same  format.   Tables  X  and  XI 
present  the  MAD  and  MDS  for  aggregate  totals  respectively. 

From  the  data  in  these  tables  it  can  be  seen  that 
lower  MADs  and  MDSs  may  result  from  increasing  the  time 
horizon.   It  has  been  suggested  that  the  lower  deviations 
were  a  mathematical  result  of  the  two-  and  three-week  time 
horizons  averaging  the  volume  fluctuations.   The  reduction 
continued  for  each  increase  in  the  time  horizon  in  8  of  16 
PODs.   In  only  one  case  did  a  decrease  not  result  from  an 
increase  in  the  time  horizon  from  one  to  two  or  three  weeks. 
In  other  words,  in  15  of  16  cases  there  was  a  reduction  in 
MAD  when  using  a  two-  or  three-week  horizon. 

The  moving  average  method  provided  a  means  of  damp- 
ening the  initial  large  volume  fluctuations  in  the  data. 
Large  volume  fluctuations  during  the  first  quarter  of  FY 
'74  could  be  dampened  by  using  a  10-week  moving  average. 
Thus,  the  first  deviation  was  not  computed  until  the  eleventh 
week. 
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The  other  methods  that  were  applied  to  the  data  be- 
gan computing  deviations  as  early  as  the  third  week.   With 
large  volume  fluctuations  at  the  beginning  of  the  period  it 
was  possible  that  there  was  a  lag  in  the  forecast  which 
initially  resulted  in  large  deviations.   The  plots  of  the 
weekly  inputs  showed  that  8  of  16  PODs  had  large  volume 
fluctuations  at  the  beginning  of  the  period.  • 

An  overall  view  of  the  forecasting  methods  indicated 
that  no  one  particular  method  dominated  the  others  in  terms 
of  providing  the  best  forecast  on  a  POD  or  an  aggregate 
level.   The  difference  between  the  largest  and  smallest  MAD 
was  less  than  100  cubic  feet  for  one  POD  (UL7) .   For  another 
POD  (UM4)  three  different  methods  were  optimal  depending  on 
the  time  horizon.   These  and  other  factors  demonstrated  the 
random  flow  of  cargo  arriving  at  the  CSS. 

6 .   Optimal  Forecasting  Method  Selection 

The  results  of  the  various  forecasting  techniques 
applied  to  the  data  provide  a  basis  for  the  selection  of  an 
optimal  method.   Through  the  application  of  these  techniques 
the  following  alternatives  were  available  for  forecasting 
the  weekly  volume  inputs  of  the  PODs : 

(1)  a  three-week  moving  average, 

(2)  a  moving  average  with  an  optimal  time  interval 
applied  to  all  PODs. 


See  Tables  X  and  XI,  Appendix  A, 
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(3)  a  moving  average  with  the  optimal  interval  for 
each  POD  applied  to  its  weekly  inputs  (i.e.,  the  time  in- 
tervals vary  among  the  PODs) , 

(4)  an  exponentially  weighted  average  with  the  op- 
timal a-level  applied  to  all  PODs, 

(5)  an  exponentially  weighted  average  with  the  op- 
timal a-level  for  each  POD  applied  to  its  weekly  inputs 
(i.e.,  the  a-level  varies  among  the  PODs), 

(6)  an  exponentially  weighted  average  with  trend 
adjustment  with  the  optimal  a-level  applied  to  all  PODs, 

(7)  an  exponentially  weighted  average  with  trend 
adjustments  with  the  optimal  a-level  for  each  POD  applied 
to  its  weekly  inputs, 

(8)  an  adaptive  response  rate  with  the  optimal 
Y-level  applied  to  all  PODs, 

(9)  an  adaptive  response  rate  with  the  optimal 
y-level  for  each  POD  applied  to  its  weekly  inputs. 

The  selection  of  the  optimal  method  required  an 
analysis  of  the  methods  used  to  calculate  the  goodness 
criterion  and  the  application  of  Student's  t-test  in  order 
to  determine  if  the  -difference  of  the  deviations  among  the 
methods  were  significantly  different. 

The  results  of  the  one  to  ten-week  moving  average 
indicated  that  a  time  interval  greater  than  three  weeks 
generated  smaller  deviations  than  the  three-week  moving 
average  for  13  of  16  PODs.   This  ratio  held  true  for  both 
criterions  and  for  all  time  horizons.   Thus,  if  a  moving 
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average  is  to  be  applied  to  the  data,  a  time  interval  greater 
than  three  weeks  is  preferred. 

There  is  one  shortcoming  in  comparing  the  moving 
average  method  with  the  others.   The  first  deviation  for  a 
moving  average  was  not  computed  until  the  week  following 
the  time  interval  of  the  average  while  the  deviations  were 
computed  beginning  the  third  week  for  the  other  methods. 
Thus,  the  longer  time  intervals  tended  to  have  a  greater 
initial  smoothing  of  the  forecasts  while  there  was  a  lag 
in  the  smoothing  for  the  other  methods.   This  fact  should 
be  kept  in  mind  when  comparing  the  aggregate  totals  in 
Table  X,  Appendix  A.   Because  of  the  manner  in  which  the 
deviations  were  computed  no  comparisons  were  made  between 
the  optimal  moving  average  and  the  other  methods. 

The  exponentially  weighted  average  offers  an  alter- 
native to  a  moving  average.   Firstly,  it  did  generate  smaller 
MADs  and  MDSs  than  a  three-week  moving  average.   A  simple 
t-test  using  MAD  as  the  criterion  and  a  one-week  time  hori- 
zon indicated  that  alternatives  (4)  and  (5)  are  significantly 
different  from  alternative  (1)  at  the  0.05  confidence  level; 
and,  therefore,  (4)  and  (5)  are  preferred  to  (1).   At  the 
same  confidence  level,  the  t-test  indicated  that  (4)  and 
(5)  are  different;  and,  therefore,  alternative  (5)  is  pre- 
ferred to  (4).   Not  until  a  confidence  level  of  0.005  is 
used,  can  it  be  said  that  alternatives  (4)  and  (5)  are  not 
significantly  different. 

The  exponentially  weighted  average  method  with  trend 
adjustment  (alternatives  (6)  and  (7))  generated  aggregate 
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MADs  and  MDSs  that  were  larger  than  those  generated  by  the 
EWA  without  the  trend  adjustment.   Based  on  this  fact,  al- 
ternatives (4)  and  (5)  are  preferred  to  alternatives  (6) 
and  (7). 

The  exponentially  weighted  average  with  an  adaptive 
response  rate  (alternatives  (8)  and  (9))  provided  another 
alternative  to  the  moving  average  method.   For  6  of  16  PODs 
the  ARR  generated  the  smallest  MAD  for  one-  and  two-week 
time  horizons.   On  the  aggregate  level  alternatives  (8)  and 
(9)  generated  slightly  smaller  MADs  than  did  alternatives 
(4)  and  (5).   A  t-test  indicated  that  alternatives  (8)  and 
(9)  are  different;  and,  therefore,  (9)  is  preferred  to  (8). 
The  test  was  conducted  at  the  0.05  confidence  level. 

With  alternatives  (5)  and  (9)  generating  smaller 
MADs  and  MDSs  than  the  three-week  moving  average,  the  ques- 
tion was  which  alternative  is  preferred?   Again,  a  simple 
t-test  at  the  0.05  level  of  confidence  was  applied.   At  this 
level,  the  test  indicated  that  alternatives  (5)  and  (9)  are 
not  significantly  different.   With  this  fact  in  mind,  al- 
ternative (5)  would  be  preferred  to  (9)  because  of  compu- 
tational  ease. 

While  it  was  possible  to  use  a  simple  t-test  to  ar- 
rive at  a  preferred  method  of  forecasting,  the  results  of 
the  various  methods  have  indicated  the  true  randomness  of 
the  cargo  flow  into  the  Container  Stuffing  Station.   No  one 


2 
Table  XII,  Appendix  A,  shows  a  sample  of  the  t-test 

calculations . 
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particular  method  dominates  the  others  in  its  ability  to 
accurately  forecast  the  weekly  volume  inputs  for  all  of  the 
major  PODs .   The  satisfaction  of  arriving  at  a  preferred 
method  must  come  from  the  fact  that  this  method  does  result 
in  deviations  that  are  smaller  than  those  generated  by  a 
three-week  moving  average  and  is  computationally  easier  and 
less  sensitive  to  the  data  than  another  method  that  is  com- 
parable in  its  forecasting  accuracy. 

Additionally,  the  minimum  MADs  achieved  using  alter- 
native (9)  occurred  at  y-levels  ranging  from  -1.00  to  +0.49. 
Also,  the  data  for  most  PODs  possessed  several  maxima  and 
minima  of  MAD  as  ylevels  were  varied.   This  revealed  the 
highly  sensitive  nature  of  the  results  of  forecasts  to  the 
input  data  when  alternative  (9)  was  used.   Should  the  char 
acteristics  of  these  data  change  in  future  observations, 
errors  resulting  from  the  continued  use  of  previously  opti- 
mal ylevels  could  become  substantial.   This  was  contrasted 
to  the  more  narrow  range  of  a-levels  (from  0.02  to  +0.77) 
for  alternative  (5)  and  the  achievement  of  a  single  minimum 
MAD  for  most  cases. 

Computer  Output  Three-One  and  Computer  Output  Five- 
One  are  representative  of  these  relationships. 
7 .   Forecast  Variance  Estimation 

To  this  point,  this  thesis  has  been  concerned  with 
the  choice  of  an  optimal  method  of  forecasting  the  volume 
inputs  to  the  CSS,  using  mean  absolute  deviation  and  mean 
deviation  squared  as  the  criterion  of  goodness.   The  appli- 
cation of  any  method  chosen  requires  that  the  statistical 
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variance  of  the  forecast  be  known  or  estimated  such  that 
confidence  in  the  forecast  may  be  assessed.   The  knowledge 
of  the  variance  can  assist  the  booking  of  vans  by  providing 
a  means  of  assessing  the  impacts  and  costs  of  conservative 
versus  aggressive  booking  policies  and  the  probable  effects 
on  average  age  of  cargo,  cancellations,  cargo  volume  on  the 
floor,  cube  utilization  and  other  policies.   Just  as  the 
expected  volume  was  forecast  by  some  method,  so  must  the 
expected  variance  be  forecasted  or  estimated.   Two  possible 
procedures  for  arriving  at  variances  for  the  forecasts  are 
suggested. 

The  first  possible  method  of  estimating  the  variance 
and  standard  deviation  of  volume  inputs  for  a  POD  could  be 
to  project  forward  historical  information  for  one  period 
to  the  next.   If  this  method  were  chosen,  the  length  of  the 
period  would  dictate  how  much  averaging  out  of  fluctuations 
would  be  accomplished.   Conversely,  should  the  fluctuations 
be  tied  to  some  major  policy  decisions  occurring  during  a 
period,  variance  estimates  could  be  seriously  in  error  if 
the  period  chosen  was  too  long.   It  is  suggested  that  per- 
haps a  period  of  six  months  to  one  year,  providing  26  to 
52  weekly  deviation  observations.   Thus,  the  weekly  volume 
input  variance  for  a  POD  calculated  for  one  period  could  be 
applied  to  all  weekly  forecasts  for  that  POD  in  the  ensuing 
period. 

Utilization  of  a  smoothed  square  deviation  computed 
simultaneously  with  each  volume  forecast  is  another  possible 
method  of  estimating  the  variance.   As  the  variance  of  each 
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forecast  is  related  to  the  square  of  the  deviation  between 
the  observed  and  forecasted  volume  for  the  period,  this 
value  could  be  smoothed  from  one  forecasting  period  to  the 
next  using  the  simple  exponentially  weighted  average  tech- 
nique.  The  estimated  standard  deviation  for  the  next  fore- 
cast could  be  calculated  as  the  square  root  of  the  smoothed 
squared  deviation.   The  value  of  the  smoothing  constant 
(alpha)  for  this  calculation  could  be  selected  depending  on 
the  weight  to  be  placed  on  the  most  recent  data  (usually 
between  0.1  and  0.3  according  to  Brown  [Ref.  3]).   Intui- 
tively, the  selection  of  a  larger  alpha  could  be  tied  to 
more  aggressive  booking  policies  by  trying  to  detect  recent 
trends  while  smaller  values  would  provide  more  conservative 
long-run  averages.   In  any  case,  this  technique  provides 
continual  updating  of  variance  estimates. 

Neither  of  these  methods  was  tested  persuant  to 
this  thesis  due  to  the  limitation  of  having  only  one  year's 
data  and  the  lack  of  any  apparent  measure  of  goodness  that 
would  indicate  the  relative  optimality  of  one  method  over 
the  other. 

Whether  or  not  either  of  these  methods  is  used, 
some  estimates  of  the  variance  associated  with  a  forecast 
prediction  is  imperative  to  establish  consistency  between 
the  booking  policy  and  other  policies  concerned  with  the 
movement  through  and  storage  within  the  Container  Stuffing 
Station . 


74 


D.   CRITIQUE  OF  PROCEDURES  AND  DATA  ANALYSIS 

The  analysis  leading  up  to  and  including  the  forecasting 
techniques  have  been  presented  in  detail  in  the  preceding 
sections  of  this  chapter.   This  section  will  evaluate  the 
procedures  used  to  make  the  analysis. 
1 .   Distribution  Audit 

The  distribution  audit  was  a  necessary  first  step 
in  the  data  analysis.   Although  there  was  a  multitude  of 
information  available  concerning  the  cargo  arriving  at  MOTBA, 
only  the  number  of  shipments,  their  volume  and  their  weight 
was  of  concern  to  this  work.   The  audit  took  the  raw  data 
and  placed  it  in  a  meaningful  format. 

The  distribution  audit  helped  to  establish  the  fact 
that  the  Vietnam  conflict  generated  a  tremendous  cargo  flow. 
Therefore,  it  was  decided  not  to  use  the  six  months  of  data 
from  FY  '73.   This  left  52  weeks  of  data  for  FY  '74,  pro- 
viding a  reasonable  number  of  observations.   Another  alter- 
native would  have  been  to  exclude  the  data  for  the  first 
two  months  of  FY  '74,  because  this  data  had  large  volume 
fluctuations  on  a  scale  similar  to  FY  '73  data.   The  ex- 
clusion of  this  data  may  have  resulted  in  the  generation  of 
lower  deviations  for  any  or  all  of  the  forecasting  techni- 
ques that  were  applied. 

Since  there  was  a  requirement  to  book  container 
space  aboard  vessels  three  weeks  in  advance,  the  data  was 
grouped  in  weekly  intervals  resulting  in  52  observations 
for  the  period.   Although  biweekly  data  could  have  been  used, 
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it  was  more  convenient  to  use  weekly  intervals  for  the  data 
analysis.   Intervals  greater  than  two  weeks  would  have  re- 
sulted in  too  few  observations. 

The  distribution  audit  also  identified  the  number 
of  shipments  classified  as  release  unit  because  of  weight. 
These  shipments  were  considered  part  of  the  random  volume 
flow  based  on  the  fact  that  they  comprised  a  very  small 
percentage  of  the  total  number  of  shipments. 

In  short,  the  distribution  audit  provided  a  start- 
ing point  for  the  procedures  that  followed.  It  simplified 
the  data  by  placing  it  in  meaningful  perspective. 

2 .  Serial  Correlation  Test 

The  data  was  analyzed  to  reveal  the  presence  of 
serial  correlation.   The  Durbin-Watson  d-Statistic  Test  is 
commonly  used  for  the  purpose.   The  calculations  for  the 
test  were  easily  accomplished  by  computer  programming.   The 
data  was  tested  using  weekly  and  biweekly  intervals.   Tri- 
weekly intervals  were  not  used  since  they  would  have  re- 
sulted in  only  17  observations.   The  results  of  the  test 
justified  the  application  of  the  forecasting  methods  pre- 
viously discussed. 

3 .  Forecasting  Methods 

a.   Moving  Average  Method 

The  moving  average  was  the  most  obvious  method 
to  initially  apply  to  the  data.   It  provided  a  simple  means 
of  forecasting  future  volume  inputs.   Less  weight  was  given 
to  the  weekly  input  as  the  time  interval  for  the  moving 
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average  was  increased.   Time  intervals  from  one  to  ten  weeks 
were  evaluated.   To  carry  the  evaluation  beyond  ten  weeks 
would  have  been  superfluous  since  time  intervals  beyond  ten 
weeks  result  in  minimal  changes  in  the  forecasts  for  each 
new  input  of  data.   As  the  moving  average  intervals  become 
large,  the  process  approaches  the  use  of  the  mean  of  the 
volume  inputs  over  a  certain  time  period  (for  example,  a 
quarter)  as  a  forecast  for  the  next  similar  period.   The 
calculations  indicated  that  increased  smoothing  of  the  data 
was  accomplished  by  increasing  the  time  interval  of  the 
moving  average  beyond  three  weeks. 

b.   Exponentially  Weighted  Average  Methods 

While  previous  discussions  may  have  indicated 
that  there  were  four  forecasting  techniques  applied  to  the 
data,  in  reality  there  were  only  two  --  the  moving  average 
and  the  EWA  plus  two  variations. 

The  EWA  method,  unlike  the  moving  average,  con- 
siders all  past  data  to  some  extent.   Through  the  applica- 
tion of  the  smoothing  constant,  the  emphasis  placed  on  the 
most  recent  data  can  be  varied. 

One  variation  of  the  EWA  was  to  adjust  for  trend 
effects  in  the  data.   The  apparent  trend  was  then  smoothed 
and  combined  with  a  smoothed  forecasted  average  to  make  a 
forecast  for  the  following  week. 

The  second  variation  was  the  application  of  an 
adaptive  response  rate  to  the  EWA.   This  particular  ARR  used 
a  smoothed  error  term  to  calculate  a  tracking  signal.   Other 
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methods  for  calculating  a  tracking  signal  may  be  available; 
but  they  were  not  investigated. 

When  the  EWA  methods  were  applied  to  the  data, 
deviations  were  calculated  beginning  with  the  third  week. 
For  the  moving  average  method,  deviations  were  not  computed 
until  the  week  following  the  length  of  the  time  interval. 
In  order  to  evaluate  the  methods  it  may  have  been  more  ap- 
propriate to  begin  calculating  the  deviations  for  the  EWA 
methods  in  the  eighth,  ninth  or  tenth  week. 
4.   Data 

The  time  period  covered  by  the  data  was  another 
limitation  in  the  application  of  the  forecasting  techniques. 
In  general,  time-series  analysis  covers  a  three  to  four  year 
period  with  data  being  grouped  by  months  or  even  quarters. 
Conditions  for  this  analysis  were  far  from  being  optimal. 

The  fact  that  there  was  only  a  year's  data  also  pre- 
cluded investigation  of  seasonal  or  cyclic  effects.   Again, 
a  three  to  four  year  period  is  usually  required  before 
smoothing  techniques  can  be  applied  to  the  data. 

In  summary,  the  breakdown  of  the  yearly  data  into 
weekly  intervals  constituted  the  most  logical  approach  to 
the  problem.   While  a  year  of  data  did  not  represent  a  long 
period  of  time,  it  did  represent  the  best  that  was  available 
at  the  time. 
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V.   CONCLUSION 

This  thesis  concerned  itself  with  the  forecast  of  the 
volume  of  cargo  flowing  into  the  Container  Stuffing  Station 
at  MOTBA  destined  for  a  particular  port  of  debarkation. 
The  data  used  in  this  work  was  taken  from  actual  operations 
at  the  CSS  during  FY  '74.   Although  every  system  is  capable 
of  generating  errors,  it  was  assumed  that  any  errors  were 
minimal  and  that  the  data  was  correct. 

Prior  to  utilizing  any  forecasting  techniques,  the  data 
was  analyzed  in  three  phases:   a  distribution  audit;  the 
grouping  of  data  into  weekly  intervals  by  POD;  and  a  test 
for  serial  correlation. 

This  thesis  made  use  of  the  distribution  audit  of 
McCarthy  and  Carter.   The  purpose  of  the  audit  was  to  sim- 
plify the  large  volume  of  raw  data  generated  by  the  opera- 
tions at  MOTBA.   The  magnitude  of  the  operations  at  the 
CSS  is  illustrated  by  the  following  information.   During 
FY  '74  there  were  a  total  of  124,592  shipments  weighing 
over  176  million  pounds  and  with  a  total  volume  of  10,070,724 
cubic  feet.   These  shipments  were  destined  for  3,356  consig- 
nees in  157  overseas  ports.   Sixteen  major  ports  accounted 
for  931  of  the  shipments,  93%  of  the  volume  and  91%  of  the 
weight.   These  16  PODs  then  became  the  focal  point  for 
further  analysis. 

The  grouping  of  the  data  into  weekly  intervals  by  POD 
indicated  the  large  volume  fluctuations  that  occurred  for 
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most  of  the  16  PODs  under  consideration.   Plots  of  the 
weekly  volume  inputs  revealed  the  lack  of  any  set  pattern 
of  cargo  flow  from  POD  to  POD.   Before  making  any  decisions 
as  to  which  forecasting  techniques  to  apply,  it  was  neces- 
sary to  determine  whether  the  data  was  random  in  nature  or 
correlated  from  one  observation  to  the  next;  that  is, 
serially  correlated.   The  Durbin-Watson  d-Statistic  Test 
satisifed  this  requirement.   A  computer  program  was  written 
to  compute  the  d-statistic  and  to  compare  it  with  critical 
values.   Using  weekly  volume  inputs  the  test  indicated  that 
9  of  16  PODs  had  serial  correlation  in  the  data.   Another 
test  was  performed  by  grouping  the  data  into  biweekly  in- 
tervals.  Similar  results  were  obtained.   The  preceding 
results  were  then  tested  to  determine  if  there  was  serial 
correlation  present  in  the  data.   The  test  was  affirmative. 

The  results  of  the  serial  correlation  test  supported 
the  use  of  forecasting  methods  other  than  ordinary  least- 
squares  trending.   The  simplest  and  most  logical  forecast- 
ing methods  to  apply  to  the  data  were  the  moving  average 
method  and  the  exponentially  weighted  average  method.   Two 
variations  were  applied  to  the  exponentially  weighted  aver- 
age resulting  in  a  total  of  four  methods. 

In  order  to  evaluate  the  forecasting  methods  on  an 
equal  basis,  the  mean  absolute  deviation  and  the  mean  de- 
viation squared  were  used,  the  deviation  being  the  differ- 
ence between  the  forecast  and  the  actual  input.   Four 
computer  programs  were  written  to  calculate  the  MADs  and 
MDSs  for  the  various  forecasting  methods. 
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Moving  averages  from  one  to  ten  weeks  were  applied  to 
the  data.   The  MAD  and  MDS  for  each  time  interval  were 
calculated  by  computer  program.   The  EWA  method  used  a 
smoothing  constant  to  generate  the  forecast.   Various  levels 
of  the  constant  were  applied  to  the  data.   A  variation  of 
the  EWA  calculated  the  apparent  trend  of  the  data.   This 
apparent  trend  and  the  forecasted  average  were  then  smoothed 
and  added  to  get  the  forecast  for  the  following  week.   An- 
other variation  to  the  EWA  method  was  the  application  of  an 
adaptive  response  rate.   The  ARR  uses  a  tracking  signal 
which  is  computed  by  taking  the  sum  of  the  differences  be- 
tween the  forecast  and  the  actual  input  volume  for  each 
week  and  dividing  this  total  into  the  sum  of  the  absolute 
values  of  those  differences.   As  each  error  becomes  avail- 
able it  is  smoothed  along  with  the  old  error.   The  smooth- 
ing constant  to  be  used  in  forecasting  is  then  set  equal  to 
the  absolute  value  of  the  tracking  signal  and  then  the  EWA 
is  applied  as  before. 

The  deviation  statistics  resulting  from  the  application 
of  the  above  forecasting  methods  to  the  actual  data  with  a 
one-week  time  horizon  were  tested  using  the  Student  t-dis- 
tribution.   This  test  revealed  that  the  EWA  and  ARR  methods 
using  the  individual  POD's  optimal  smoothing  constants  were 
statistically  the  same  and  better  than  all  other  methods. 
The  EWA  technique  was  then  chosen  as  preferred  based  on 
its  computational  ease  relative  to  ARR. 

The  selection  of  the  EWA  method  offered  some  advantages 
over  the  three-week  moving  average  being  used  at  the  CFD. 
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First,  EWA  was  computationally  no  more  difficult.   It  de- 
rived its  advantage  by  responding  to  changes  in  volume  more 
uniformly  than  the  abrupt  alterations  that  the  three-week 
moving  average  caused  in  the  forecasts.   Also,  the  accumula- 
tion of  a  historical  data  base  in  the  smoothed  forecast 
average  increased  the  accuracy  of  the  EWA  method  as  the  time 
horizon  increased.   Decreases  of  4.91,  8.0%  and  10.01  in 
mean  absolute  deviations  for  the  one-  ,  two- ,  and  three-week 
time  horizons  using  the  EWA  method  evidenced  the  increase 
in  accuracy.   For  this  application  and  data,  there  can  be 
little  doubt  that  the  use  of  the  EWA  provided  improvement 
in  the  forecasting  techniques. 

The  other  forecasting  methods  retained  their  relative 
merits.   Greater  accuracy  was  achieved  using  moving  averages 
by  increasing  the  lengths  of  the  averaging  interval  at  the 
expense  of  increasing  the  computational  difficulty.   The 
absence  of  any  apparent  trend  in  the  data  resulted  in  lack 
of  any  additional  benefits  from  trend  adjustments.   While 
the  adaptive  response  rate  technique  showed  no  real  advan- 
tage over  the  EWA  method  for  this  data,  its  ability  to 
respond  rapidly  to  major  policy  changes  makes  it  a  viable 
alternative  at  the  expense  of  increased  computational  dif- 
ficulty.  As  more  data  becomes  available,  all  of  these 
methods  deserve  additional  consideration  and  thoughtful 
attention  as  the  results  of  this  work  are  based  on  only  a 
year's  data. 

Although  significant  increases  in  forecasting  have  been 
achieved,  there  remains  in  the  data  large  amounts  of  randomness 
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causing  significant  errors.   Its  random  nature,  without 
trends  or  significant  shifts,  have  led  to  the  evolution  of 
the  exponentially  weighted  average  method  as  the  preferred 
choice.   As  more  data  is  accumulated  the  other  methods  may 
be  reexamined.   The  optimal  smoothing  constants  should  be 
updated  periodically  to  account  for  changes  in  the  data  and 
promote  long  run  efficiency  of  the  method.   However,  the 
most  important  safeguard  against  costly  mistakes  due  to 
random  effects  in  the  input  volume  is  careful  managerial 
scrutiny  of  each  forecast  generated  and  used  to  make  a 
booking.   Judgment  on  the  reasonableness  of  a  prediction 
makes  the  manager  the  final  and  ultimate  smoothing  factor 
in  the  utilization  of  the  forecast. 
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APPENDIX 

A 

WEEK 

SHIPMENTS 

VOLUME 

WEIGHT 

27 

3,208 

220,991 

3,925,310 

28 

5, 

,091 

413, 

,602 

7 

,164. 

,543 

29 

4, 

598 

336. 

,189 

6, 

,062 

,987 

30 

4, 

884 

368, 

,350 

6. 

,038 

,153 

31 

5 

,707 

344. 

,109 

6, 

312 

,502 

32 

4. 

,727 

372. 

,533 

6 

,281, 

,837 

33 

4, 

726 

358, 

880 

6. 

,640. 

,755 

34 

4. 

,214 

303. 

,929 

5 

,393 

,130 

35 

5. 

,681 

378. 

598 

7. 

,275 

,388 

36 

5 

,031 

324, 

176 

6 

,206 

,201 

37 

5, 

045 

348. 

497 

6. 

,760 

,475 

38 

4 

,638 

301, 

952 

5 

,808 

,959 

39 

4, 

384 

303. 

536 

5. 

,719 

,112 

40 

4. 

,047 

300. 

089 

5 

,138 

,948 

41 

4, 

646 

335, 

702 

6. 

322 

,970 

42 

4, 

753 

343. 

,370 

6 

,425 

,117 

43 

4. 

,317 

314. 

323 

5 

,969 

,897 

44 

3. 

,935 

281, 

744 

5 

,190 

,323 

45 

3. 

906 

253. 

313 

4 

,840 

,206 

46 

3 

,318 

248 

,060 

4 

,473 

,395 

47 

3 

660 

288 

,476 

5 

085 

827 

48 

2. 

957 

231, 

423 

4 

,205 

,199 

49 

3 

,498 

320 

,646 

5 

,902 

,697 

50 

3 

,338 

285 

413 

1 

,150 

,959 

51 

3 

,675 

387 

,533 

5. 

,980 

,555 

52 

3 

,579 

433 

,497 

6 

,248 

,080 

Note:   (1)   Volume  is  in  cubic  feet 
(2)   Weight  is  in  pounds 

Table  I.   Gross  Weekly  Inputs.   Weeks  27  through  52,  FY  '73 
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WEEK 

SHIPMENTS 
2,444 

VOLUME 
261,448 

WEIGHT 

1 

3,921,219 

2 

2,797 

287, 

,444 

4, 

,527, 

403 

3 

3,141 

278, 

991 

4. 

,997. 

,954 

4 

2,758 

232. 

,296 

3 

,571, 

542 

5 

3,414 

284, 

814 

5j 

151, 

369 

6 

4,044 

308, 

446 

5, 

,844 

,244 

7 

965 

88. 

,024 

1. 

,484, 

204 

8 

2,905 

208. 

,179 

3. 

,908 

,164 

9 

2,641 

196 

,694 

3. 

,625. 

615 

10 

1,994 

164 

,170 

3 

,017 

,826 

11 

2,425 

206. 

,833 

3 

,837, 

,876 

12 

2,787 

215. 

422 

4 

,308. 

528 

13 

2,485 

199. 

,702 

4 

,533 

,820 

14 

2,541 

192. 

863 

3 

,600 

,097 

15 

1,784 

135, 

561 

2 

,378. 

,525 

16 

2,653 

194, 

927 

3 

,742 

,275 

17 

1,983 

150. 

,028 

2 

,857. 

,312 

18 

1,387 

HI. 

,154 

2 

,199. 

,669 

19 

2,340 

208. 

,560 

3. 

,761 

,167 

20 

2,411 

184 

,520 

3. 

,200 

,648 

21 

1,361 

105, 

629 

1 

,965. 

,960 

22 

2,414 

149. 

722 

3. 

,663 

,732 

23 

2,927 

203 

,497 

3 

,864. 

,209 

24 

2,754 

199, 

888 

3 

,695 

,825 

25 

2,757 

189 

,969 

3 

,657 

,957 

26 

1,523 

130. 

513 

2 

,592 

,805 

Note:   (1)   Volume  is  in  cubic  feet 
(2)   Weight  is  in  pounds 


Table  II.   Gross  Weekly  Inputs.   Weeks  1  through  26,  FY  » 74 
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WEEK 

SHIPMENTS 

VOLUME 

WEIGHT 

27 

1,530 

135,843 

2,186,285 

28 

2,956 

232. 

,446 

3, 

841,574 

29 

3,058 

260, 

,123 

4. 

,468,724 

30 

2,338 

179, 

400 

3, 

078,738 

31 

2,580 

241 

,953 

3. 

793,319 

32 

2,814 

242 

,661 

4. 

,034,761 

33 

2,667 

211 

,740 

3. 

,470,228 

34 

1,585 

139, 

665 

2. 

,459,631 

35 

3,007 

267. 

,702 

4, 

907,835 

36 

2,083 

164. 

,209 

2. 

,725,845 

37 

1,994 

143. 

,360 

2 

,385,097 

38 

2,198 

145 

,621 

2. 

,851,101 

39 

2,557 

172. 

,236 

3. 

,151,837 

40 

2,247 

170. 

,630 

2 

,921,166 

41 

2,060 

201 

,845 

3 

,388,849 

42 

2,145 

195. 

,070 

3 

,211,454 

43 

2,041 

191 

,623 

3 

,135,205 

44 

2,704 

198 

,084 

3 

,052.520 

45 

2,381 

204 

,956 

3 

,269,329 

46 

2,189 

193 

,582 

3 

,264,435 

47 

2,739 

177 

,121 

3 

,355,770 

48 

1,395 

96 

,791 

1 

,628,575 

49 

2,227 

185 

,867 

2 

,771,581 

50 

2,522 

194 

,153 

2 

,716,456 

51 

2,812 

229 

,928 

3 

,602,602 

52 

2,228 

188 

,520 

2 

,631,603 

Note:   (1)   Volume  is  in  cubic  feet 
(2)   Weight  is  in  pounds 


Table  II.   Gross  Weekly  Inputs 
(Continued) 


Weeks  27  through  52,  FY  '74 
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Period  analyzed  -  FY  '74 

Shipments 

Volume  of  cargo 
Weight  of  cargo 
Number  of  PODs 

Number  of  consignees 


365  days 

124,592  shipments 

10,070,724  cubic  feet 

176,213,333  pounds 
157 

3,356 


Table  III.   Data  Summary  for  the  Distribution  Audit. 
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SHIPMENTS 

Total  shipments  of  16  PODs 
Total  shipments  of  all  PODs 
Percentage  of  total  shipments 


115,732  shipments 
124,592  shipments 
92.891 


VOLUME 

Total  volume  of  16  PODs 

Total  volume  of  all  PODs 

Average  volume/shipments  (16  PODs) 

Percentage  of  total  volume  (16  PODs) 


9,370,693  cubic  feet 

10,070,724  cubic  feet 

87  cubic  feet 


93.05 


WEIGHT 

Total  weight  of  16  PODs 

Total  weight  of  all  PODs 

Average  weight/shipments  (16  PODs) 

Percentage  of  total  weights  (16  PODs) 


160,575,726  pounds 

176,213,333  pounds 

1,387.47  pounds 
91% 


Table  IV.   Tabular  Data  for  16  Major  PODs,  FY  '74. 


POD 

SHIPMENTS 

VOLUME 

WEIGHT 

RA3 

20,125 

1,474,747 

28,284,704 

RGU 

6,041 

491,465 

11,486,615 

RG1 

1,850 

170,766 

3,774,348 

RJ1 

2,848 

170,635 

3,323,333 

RJ3 

2,210 

106,150 

1,802,072 

SA3 

16,837 

1,377,838 

19,920,384 

TA1 

2,502 

120,073 

2,132,222 

TA2 

6,793 

502,989 

7,269,693 

UB1 

12,339 

1,081,519 

19,311,360 

UC2 

6,593 

976,869 

19,325,856 

UD6 

8,691 

451,466 

9,062,837 

UL7 

3,183 

. 233,232 

4,445,639 

UM1 

6,619 

434,615 

7,081,348 

UM4 

3,769 

276,100 

3,740,452 

XE2 

14,008 

1,315,539 

16,679,435 

XJ1 

1,324 

186,690 

2,935,428 

Note:   (1)   Volume  is  in  cubic  feet 
(2)   Weight  is  in  pounds 


Table  V.   POD  Data  -  16  PODs ,  FY  '74 
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NUMBER  OF  PERCENTAGE  OF 

POD  RU  SHIPMENTS              RU  SHIPMENTS 

RA3  384  1.90 

RGU  156  2.58 

RG1  37  2.00 

RJ1  41  1.43 

RJ3  15  0.68 

SA3  183  1.00 

TA1  15  0.59 

TA2  86  1.26 

UB1  242  1.96 

UC2  327  4.90 

UD6  80  0.92 

UL7  43  1.35 

UM1  80  0.92 

UM4  20  0.53 

XE2  136  0.97 

XJ1                    5_5  4.15 

TOTAL  1,900 

The  percentage  of  RU  shipments  to  the  total  number  of 
shipments  for  the  16  major  PODs  is  1.64%. 


Table  VI.   Release  Unit  Shipments 
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WEEK 
POD  INT 


RA3 


WEEK  NO.  d 
START  OBS 


52 
26 
25 


STATISTIC 

MEANING 


INCN 
INCN 
POSI 


WEEK 
POD  INT 


UB1 


WEEK   NO. 
START  OBS 


52 
26 

25 


d-STATISTIC 

MEANING 

NONE 
NONE 
NONE 


RGU 


RG1 


RJ1 


RJ3 


SA3 


TA1 


TA2 


52 
26 

25 

52 
26 

25 

52 
26 
25 

52 
26 
25 

52 
26 

25 

52 
26 

25 

52 
26 

25 


POSI 

INCN 
INCN 

POSI 
POSI 
POSI 

POSI 
NONE 
NONE 

NONE 
NONE 
NONE 

POSI 
INCN 
POSI 

INCN 
POSI 
INCN 

INCN 
INCN 
INCN 


UC2 


UD6 


UL7 


UM1 


UM4 


XE2 


XJ1 


52 
26 
25 

52 
26 

25 

52 
26 

25 

52 
26 

25 

52 
26 
25 

52 
26 
25 

52 
26 
25 


NONE 
INCN 
INCN 

INCN 
NONE 
NONE 

POSI 
POSI 
POSI 

POSI 
NONE 
POSI 

POSI 

NONE 
POSI 

POSI 
POSI 
POSI 

POSI 
POSI 
INCN 


Note:   (1)  POSI  =  the  test  indicated  POSITIVE  serial  correlation 

(2)  INCN  =  the  test  was  INCONCLUSIVE 

(3)  NONE  =  the  test  indicated  no  serial  correlation 


Table  VII.   Serial  Correlation  Test  Results 
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Table  VIII.   MADs  for  16  PODs . 

One-week  time  horizon 


POD 

TWMA 

OMA 

EWA 

EWA-TA 

EWA 

-ARR 

t 

a 

a 

1 

RA3 

6797* 

3 

6797* 

.36 

7021 

.11 

7186 

-  .78 

7074 

RGU 

3910 

6 

3733 

.27 

3599* 

.13 

3717 

.40 

3633 

RG1 

1754 

10 

1467* 

.60 

1709 

.34 

1787 

.34 

1611 

RJ1 

1793 

9 

1462* 

.35 

1734 

.14 

1770 

.01 

1604 

RJ3 

1130 

8 

1064 

.21 

1083 

.11 

1066 

-  .01 

1030* 

SA3 

7862 

3 

7862 

.04 

7458* 

.02 

7461 

-  .05 

7974 

TA1 

996 

5 

811 

.22 

935 

.10 

961 

.13 

395 

TA2 

3654 

10 

3498 

.40 

3505 

.18 

3469* 

-  .02 

3484 

UB1 

5321 

7 

5040* 

.24 

5112 

.10 

5174 

-  .85 

5216 

UC2 

5236 

7 

4583* 

.10 

4940 

.04 

4839 

-  .17 

4794 

UD6 

4205 

7 

3718 

.28 

3845 

.11 

3893 

-  .07 

3372* 

UL7 

1383 

5 

1296 

.20 

1290 

.  12 

1358 

.38 

1287* 

UM1 

3637 

10 

3341 

.54 

3300 

.10 

3311 

-  .10 

3161* 

UM4 

2521 

8 

2366 

.48 

2380 

.21 

2490 

.11 

2355* 

XE2 

7770 

10 

7391 

.47 

7615 

.26 

7787 

.07 

7130* 

XJ1 

2743 

2 

2564 

.43 

2434* 

.21 

2657 

-1.00 

2475 

Note:   (1)  MAD  =  Mean  Absolute  Deviation 

(2)  POD  =  Port  of  Debarkation 

(3)  TWMA  =  Three  Week  Moving  Average 

(4)  OMA  =  Optimal  Moving  Average 

(5)  t  =  time  interval  (in  weeks)  for  OMA 

(6)  EWA  =  Exponentially  Weighted  Average 

(7)  a  =  smoothing  constant  for  EWA  and  EWA-TA 

(8)  EWA-TA  =  Exponentially  Weighted  Average  with 

Trend  Adjustment 

(9)  EWA-ARR  =  Exponentially  Weighted  Average  with 

Adaptive  Response  Rate 

(10)  y  -    smoothing  constant  for  EWA-ARR 

(11)  MADs  are  expressed  in  cubic  feet 

*Indicates  smallest  MAD  for  POD 
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Table 

VIII. 

MADs  for  16 

PODs 

(Continued) . 

Two- 

•week  time 

horizon 

POD 

TWMA 

t 

OMA 

EWA 
a  

a 

EWA-TA 

EWA-ARR 
I  

RA3 

6204 

3 

6204 

.32 

5955* 

.17 

6293   - 

.89 

6095 

RGU 

3621 

3 

3621 

.09 

32.84 

*.22 

3489 

.45 

3339 

RG1 

1601 

10 

1362 

.77 

1572 

.37 

1690 

.23 

1311* 

RJ1 

1350 

8 

1194* 

.37 

1420 

.17 

1449 

.01 

1292 

RJ3 

850 

3 

850 

.28 

804 

.11 

801 

.01 

750* 

SA3 

6639 

3 

6639 

.02 

6250 

.01 

6247*  - 

.05 

6927 

TA1 

798 

6 

595* 

.22 

701 

.12 

729 

.11 

694 

TA2 

3194 

10 

3056 

.39 

3072 

.17 

3125   - 

.03 

2991* 

UB1 

4655 

7 

4250 

.15 

4246 

.06 

4179*  - 

.01 

4417 

UC2 

4001 

7 

3150* 

.15 

3399 

.06 

3278   - 

.10 

3423 

UD6 

3688 

7 

3135 

.27 

5255 

.10 

3268   - 

.08 

2938* 

UL7 

1053 

4 

1005* 

.28 

1013 

.17 

1082   - 

.44 

1033 

UM1 

3146 

10 

2641 

.19 

2816 

.09 

2893 

.01 

2615 

UM4 

2343 

9 

2153  - 

.02 

1878 

-.01 

1847* 

.10 

2177 

XE2 

6984 

3 

6984 

.65 

6697 

.31 

6941   - 

.01 

6269* 

XJ1 

2771 

10 

2327* 

.23 

2435 

.10 

2592   - 

.75 

2452 

*Indicates  smallest  MAD  for  POD 
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Table  VIII.   MADs  for  16  PODs  (Continued). 

Three-week  time  horizon 


POD 

TWMA 

OMA 

EWA 

EWA-TA 

EWA-ARR 

t 

a 

a 

I 

RA3 

6620 

8 

6197 

.26 

5910 

.10 

6069 

-  .74 

6351 

RGU 

3477 

5 

3410 

.13 

3120* 

.04 

3382 

-  .47 

3214 

RG1 

1652 

10 

1363 

.75 

1498 

.37 

1529 

.25 

1245* 

RJ1 

1188 

7 

1036* 

.43 

1253 

.25 

1296 

.01 

1086 

RJ3 

863 

10 

824 

.31 

789 

.09 

784 

.01 

713* 

SA3 

6451 

7 

6388 

.34 

5848 

.01 

5238* 

-1.00 

6234 

TA1 

721 

5 

545* 

.25 

626 

.13 

653 

.14 

561 

TA2 

3149 

10 

2660* 

.31 

3030 

.17 

2990 

-  .03 

2789 

UB1 

3627 

7 

3252 

.18 

3225 

.07 

3187* 

-  .01 

3471 

UC2 

3390 

10 

2548* 

.19 

2830 

.06 

2781 

-1.00 

2652 

UD6 

3326 

10 

2981 

.29 

2988 

.11 

3147 

-  .07 

2410* 

UL7 

972 

5 

937* 

.29 

946 

.16 

995 

-  .46 

978 

UM1 

3215 

10 

2355 

.19 

2675 

.09 

2681 

-  .12 

2331* 

UM4 

2327 

9 

1970  - 

-.01 

1709* 

.01 

2051 

-  .01 

2127 

XE2 

6831 

10 

6717 

.55 

6778 

.30 

6960 

-  .01 

5971* 

XJ1 

2787 

10 

2049* 

.17 

2308 

.09 

2422 

.75 

2386 

*Indicates  smallest  MAD  for  POD 
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Table  IX.   MDSs  for  16  PODs 
One-week  time  horizon 


POD 

RA3 
RGU 
RG1 
RJ1 
RJ3 
SA3 
TA1 
TA2 
UB1 
UC2 
UD6 
UL7 
UM1 
UM4 
XE2 
XJ1 

Note: 


TWMA 


92,564  7 

27,737  4 

10,349  10 

5,741  8 

2,825  7 
96,136*  3 

1.659  6 

22,695  10 

51,088  7 

41,418  9 

34,556  7 

3,453  5 

25,750  10 

11,870  8 

102,515  10 

15,387  2 


OMA 

90,250* 

27,546 
5,769* 
3,691* 
2,591 
96,136* 
1,184* 
19,113* 
41,265* 
30,695* 
30,086 
3,195 
23,983 

9,145* 
84,346* 
13,530 


EWA 


EWA-TA 


a 


30  93 

32  24 
60  10 

33  5 
12  2 
38  98 
21  1 
28  20 
14  42 
14  35 
25  29 
28  3 
24  22 
28  9 
52  103 
48  12 


854 

977* 

736 

337 

455* 

728 

545 

456 

697 

775 

872 

183 

108 

936 

848 

845 


a 


06  98,367 

14  26,690 

30  10,117 

16  5,511 

05  2,481 

04  104,177 
12  1,611 

12  21,253 

06  43,547 

05  34,811 

13  31,267 
16  3,267 
11  23,010 

07  10,610 
25  106,797 
16  14,115 


EWA-ARR 
X 

-  .82  98,528 

-  .59  25,060 
.25  8,805 
.02  5,474 

-  .55  2,455 
-1.00  104,207 

.11  1,614 

-  .03  20,549 

-  .78  44,880 
-1.00  35,382 

-  .07  27,730* 

-  .44  3,176* 

-  .14  21,463* 

-  .86  10,500 
.05  90,932 

-  .13  12,490* 


(1)  MDS  =  Mean  Deviation  Squared 

(2)  POD  =  Port  of  Debarkation 

(3)  TWMA  =  Three  Week  Moving  Average 

(4)  OMA  =  Optimal  Moving  Average 

(5)  t  =  time  interval  (in  weeks)  for  OMA 

(6)  EWA  =  Exponentially  Weighted  Average 

(7)  a  =  smoothing  constant  for  EWA  and  EWA-TA 

(8)  EWA-TA  =  Exponentially  Weighted  Average  with 

Trend  Adjustment 

(9)  EWA-ARR  =  Exponentially  Weighted  Average  with 

Adaptive  Response  Rate 

(10)  y  =  smoothing  constant  for  EWA-ARR 

(11)  MDSs  are  in  1,000*5  and  expressed  in  (cubic  feet) 

indicates  smallest  MDS  for  POD 
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Table  IX.   MDSs  for  16  PODs 
(continued) 


Two-week  time  horizon 


POD   TWMA 


OMA 


EWA 


EWA-TA 


a 


a 


RA3  70,264  7  69,626 

RGU  21,534  5  20,533 

RG1  9,064  10  5,293* 

RJ1  3,732  8  2,413* 

RJ3  1,659  10  1,388 

SA3  76,062  3  76,062 

TA1  1,060  6  749* 

TA2  16,801  10  12,717* 

UB1  30.872  7  24.900* 

UC2  25,069  10  14,927* 

UD6  24,416  10  20,683 

UL7  2,084  4  1,880* 

UM1  19,208  10  15,958 

UM4  10,062  9  6,781* 

XE2  81,399  2  80,557 

XJ1  13,673  10  10,908* 


.18  62,116* 

.27  18,672* 

.69  8,385 

.32  4,019 

.13  1,327* 

.07  73,856* 

.22  953 

.26  14,376 

.14  25,001 

.15  18,843 

.23  21,071 

.30  1,952 

.20  16,063 

.15  7,909 

.54  79,063 

.31  11,087 


05  64,218 

12  20,411 
36  7,121 

16  4,165 

06  1,360 
03  73,992 

13  1,018 

11  15,107 
06  25,749 
06  17,843 

12  22,598 

17  2,021 
10  16,948 
01  8,057 
26  81,152 

13  12,110 


EWA-ARR 
I 

-  .82  69,925 

-  .53  18,829 
.22  5,588 
.01  3,965 

-  .55  1,334 
-1.00  80,384 

.11  948 

-  .03  13,101 

-  .01  27,513 
-1.00  17,496 

-  .07  17,668* 

-  .47  1,977 

-  .01  14,948* 

-  .01  8,216 

-  .01  60,822* 

-  .14  11,067 


*Indicates  smallest  MDS  for  POD 
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Table  IX.   MDSs  for  16  PODs 
(continued) 

Three-week  time  horizon 


POD 

TWMA 

t 

OMA 

a 

EWA 

a 

EWA-TA 

I 

EWA-ARR 

RA3 

69,762 

8 

64,618 

.10 

54,579* 

.04 

55,583 

-  .78 

67,596 

RGU 

19,006 

4 

18,391 

.26 

16,418* 

.11 

18,245 

-  .51 

16,570 

RG1 

10,187 

10 

5,461* 

.75 

8,571 

.38 

6,724 

.24 

5,822 

RJ1 

2,737 

10 

1,871* 

.33 

3,234 

.18 

3,354 

-  .01 

2,968 

RJ3 

1,347 

9 

1,035 

.14 

967* 

.06 

1,002 

-  .55 

969 

SA3 

74,786 

3 

74,786 

.05 

61,446 

.02 

61,394* 

-1.00 

73,494 

TA1 

806 

10 

663 

.26 

712 

.14 

779 

.11 

635* 

TA2 

15,444 

10 

9,783* 

.24 

12,773 

.10 

13,408 

-  .03 

10,545 

UB1 

18,965 

4 

16,064 

.17 

15,866* 

.07 

16,757 

-  .01 

17,605 

UC2 

19,213 

10 

10,170* 

.16 

13,467 

.06 

12,539 

-1.00 

11,723 

UD6 

20,080 

10 

15,516 

.24 

16,499 

.13 

18,160 

-  .07 

12,549* 

UL7 

1,633 

4 

1,551* 

.33 

1,570 

.18 

1,638 

.49 

1,621 

UM1 

17,862 

10 

12,662 

.19 

13,577 

.10 

14,584 

-  .12 

11,984* 

UM4 

8,787 

8 

5,197* 

.01 

6,067 

.01 

6,217 

-  .69 

7,084 

XE2 

81,188 

10 

75,150 

.53 

79,907 

.26 

80,197 

-  .01 

54,919* 

XJ1 

13,778 

10 

8,106* 

.23 

10,803 

.12 

11,790 

-  .94 

10,820 

*Indicates  smallest  MDS  for  POD 
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WEEKLY 

TIME      TWMA       OMA         EWA  EWA-TA  ARR 

HORIZON   t a a y_ 

one     545,749   7  523,120  .37  535,979  .17  566,904  -.14  543,745 

two     406,966  10  391,469  .36  388,940  .17  421,432  -.05  393,060 

three   375,569  10  341,400  .36  350,451  .17  384,810  -.05  346,852 


Note : 

(1)  MDS  =  Mean  Deviation  Squared 

(2)  OMA  =  Optimal  Moving  Average  (a  specific  optimal  interval 
is  applied  to  all  POD  data) 

(3)  EWA  =  Exponentially  Weighted  Average  (a  specific  optimal 
a-level  is  applied  to  all  POD  data 

(4)  EWA-TA  =  Exponentially  Weighted  Average  with  Trend  Ad- 
justment (a  specific  a-level  is  applied  to  all  POD  data) 

(5)  ARR  =  Adaptive  Response  Rate  (a  specific  optimal  ylevel 
is  applied  to  all  data 

(6)  t  =  optimal  time  interval  for  OMA 

(7)  a  =  smoothing  constant  for  EWA  and  EWA-TA 

(8)  y  =  smoothing  constant  for  ARR 

(9)  MDSs  are  expressed  in  (cubic  feet)2 


Table  XI.   Aggregate  MDSs  for  16  PODs . 
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In  this  work  the  Student  t-distribution  was  used  to 
resolve  whether  the  mean  of  the  MAD  differences  between  two 
forecasting  methods  was  equal  to  zero. 

The  hypotheses  were: 

HQ:  u  =  0  ,   Hx:  u  +   0. 

The  t-test  utilizes  the  following 

t  =  (X  -  u)  /n" 
s 

where  X  is  the  sample  mean,  u=0,  n=number  of  observations 

and  s=the  sample  standard  deviation. 

Calculated  values  of  t  are  compared  with  critical  values 

of  t.   If 

t  (calculated)  <  t  (critical) 

H   is  accepted 
o        r 


If 


t  (calculated)  >  t  (critical) 

H   is  not  accepted. 

o  v 


If  H   is  not  accepted,  the  two  methods  are  different  and 
o  v         ' 

the  method  that  generates  the  lowest  mean  absolute  deviation 
is  chosen  as  the  preferred. 

The  test  only  utilized  deviations  generated  by  a  one- 
week  time  horizon. 


Table  XII.   t-Test  Procedure 
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COMPUTER   OUTPUT    ONE 


TABLE    Of    STAT 

POO    INTEKVAL 

1 
2 
2 


1SIICAL    DATA    FOR 
START 


SFLfcCTcD    POOS    USING    WEEKLY     TOTAL    VOLUME     (CIJ.     FT.)     DATA 


AC  2 


AF4 


8A1 


8A3 


«JA4 


CA2 


C/.3 


CB5 


Ctl 


CE2 


CHI 


CJ1 


CKA 


OKI 


CK2 


EOl 


ED3 


E04 


ERCtPT 

SLOPE 

STu.EKR. 

R-SOUARE 

O-STAT 

08  S. 

MEANING 

0. 

125. 

10. 

0.4l4o7E 
0  .  1  7  a  1  '>  E 

0  .  8  2  7  0  1  f 

DO 
01 
01 

1  75.1 

247.5 
239.  > 

0.00131 
0.00313 
0.06791 

1  .60832 
1.  34058 
1  .'?0o06 

52 
26 
25 

MCN5 
INCN 

NONE 

2. 
5. 
7. 

0.99676E- 
0.353i.5E 
0.311 7oE 

-01 
00 

00 

2R.4 
40.2 
40.3 

0.00289 
0.00469 
0.00363 

2.0728R 
2.14859 
2.16524 

5? 
26 
25 

NON^ 
NONE 
NONE 

56. 

115. 

40. 

-0.11562E 

-O.48520E 
-0.4841  IE- 

01 

01 

-01 

133.9 
167.  7 
183. 2 

0.01717 

0.03^14 
0.00223 

1  .99740 
1  .11634 
0.98535 

52 
26 
?5 

NONE 
POSI 
POSI 

973. 
1936. 
1901. 

-0.  141 74E 
-0.4H123E 
-0.6625 IE 

01 
01 
01 

1588.2 
24  32.2 
23o0.8 
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0.0002>> 
0.05028 

1 . 79631 
1.75145 
1  .32076 

52 

26 
25 

NOME 

NONE 
NONE 

3. 

c . 
220. 

-0.576i8E 
-0.21539E 

-0. 12810E 

-01 
00 
02 

12.6 

1  7.9 

277.3 
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2.0 *8f3 
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25 

NONE 
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184. 
3/3. 
.4*1. 
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-0.131 9 7E 

01 
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32 
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0.4 
6.  c 

J. 02979 
O.O61 18 
0.11932 

2.101 97 
2. 22>12 
0.1 14* 3 

5? 
26 
23 

NONE 
NONE 
POSI 

25. 

50. 
4b. 

-0.36933E 
-0. 14968E 
-0.12378E 

00 
01 
01 

32.1 

46.6 
50.  1 

0.03002 
0.059?3 
0.03582 

1.83209 
1.64072 
1.48792 

52 
26 
25 

NONE 
NONE 
NONE 

101 


COMPUTER  OUTPUT   TWO 
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COMPUTER   OUTPUT   THREE -ONE 
Mean  Absolute   Deviation    (in   cubic    feet) 
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0.40 

59477. 

536504832. 

0.39 

59437. 

536179968. 

0.38 

59409. 

536002560. 

0.3/ 

59388. 

535979264. 

0.36 

593o7. 

536119552. 

0.35 

59402. 

536432123. 

0.34 

59424. 

536930304. 

0.33 

59461. 

537623296. 

0.32 

59505. 

538528000. 
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COMPUTER   OUTPUT    FOUR-ONE 
Mean  Absolute   Deviation    (in   cubic    feet) 
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COMPUTER   OUTPUT    FOUR-TWO 
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COMPUTER  OUTPUT  FOUR-THREE 

ALPHAS 


TABLE  CF 

AGGREGATE 

DEVIATIONS  FOR  ALL 

1.00 

118467. 

2101308416. 

0.99 

117072. 

2C49d01728. 

C.98 

115695. 

2000065024. 

0.97 

114331. 

195201740b. 

C.96 

112984. 

1905584384. 

0.95 

111655. 

1860696624. 

0.94 

110342. 

1817293056. 

0.93 

109076. 

177530342<t. 

C.92 

107855. 

1734674176. 

0.91 

106647. 

1695347968. 

C.90 

105467. 

1657269760. 

0.89 

104315. 

1620392443. 

0.68 

103180. 

1584663552. 

0.87 

102075. 

1550044928. 

0.86 

100982. 

151648V472. 

0.85 

99902. 

1463955968. 

0.84 

98834. 

1452407808. 

0.83 

S7789. 

14218C4544. 

0.82 

96763. 

1392116736. 

0.81 

95766. 

1363305216. 

0.80 

94803. 

1335342592. 

0.79 

93866. 

1308194816. 

0.78 

92932. 

1281838336. 

0.77 

92005. 

1256241920. 

0.76 

91083. 

1231379712. 

0.75 

90166. 

12C7229696. 

0.74 

89262. 

1183765248. 

0.73 

88388. 

1160964096. 

0.72 

87540. 

1138805504. 

0.71 

86729. 

1117269504. 

0.70 

85948. 

1096335872. 

0.69 

85188. 

1075986176. 

0.68 

64^53 • 

1056201984. 

0.67 

83765. 

1036966912. 

0.66 

83084. 

1018264832. 

0.65 

82419. 

1000078592. 

0.64 

81775. 

982395643. 

0.63 

81141. 

965200384. 

C.62 

80502. 

948480000. 

0.61 

75859. 

932221440. 

0.60 

79219. 

916412923. 

0.59 

78581. 

901042176. 

0.58 

77959. 

886095104. 

0.57 

77368. 

871566592. 

0.56 

76803. 

857441792. 

0.55 

76237. 

843713024. 

0.54 

75679. 

830370048. 

0.53 

75125. 

817404416. 

0.52 

74571. 

804806656. 

0.51 

74C29. 

792569088. 

0.50 

73495. 

780682496. 

0.49 

72962. 

769141248. 

0.48 

72429. 

757939456. 

0.47 

71901. 

747066368. 

0.46 

71338. 

736519163. 

0.45 

70871. 

726268896. 

0.44 

70363. 

716374784. 

0.43 

69850. 

706766336. 

0.42 

69335. 

697462016. 

0.41 

68835. 

688456448. 

0.40 

68339. 

679745536. 

0.39 

67840. 

671324928. 

0.38 

67356. 

663196672. 

0.37 

66887. 

655353856. 

0.36 

66416. 

647796736. 

0.35 

65952. 

640525312. 

0.34 

65493. 

633539584. 

0.33 

65042. 

626338784. 

0.32 

64616. 

620430080. 

COMPUTER  OUTPUT  FIVE-ONE 


Mean  Absolute    Deviation    (in   cubic    feet) 
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COMPUTER  OUTPUT    FIVE-TWO 


C9     rsuj^-rryocga^rgrn^-ogrg?"  um^^Tvr1/1  t  ro  rg  o  o  <  o  ^a  rn  o 'J*  orgcr>f\JCr'C7*fyieo  0~-4O<^Jr^pgo  r  oin^oo«X)'.7»*r<-r  ^  4"Orr»i7'— •mo'  oomfns'j 
in     r-  m^o^^a  on  gr  g>  ^O'o  jToincou^f^-co  oco.r\,"OiA^rgr\jorrtfVJf*""  oo  o<  vj^^n  g>  'n  o<r\f>ir^-u^-^Lrn/»r^.oO''Orno  or-o*-nrvj— «x)-^r^-rjrsj— <  f  co  r 

O      CONgj^POCNJ-^OON-Otor-'A  T*>f^(NJ~i-<oOO^^^'C^^(7,OOO^Hf\jrvjln^^«fi/\^  O  O^-f^iD  OCrO^OOOJ1  Jv^y*'i)ft)'0'T,O^Nr^'ni)'T) 

a.    w"»-Au~nAu~»iAiAu^gr  g^grgfgrvr.rgj'J/grgrgrgr^rrnrornrnrnrnr^g'grg^^^grvrg^grgr-Tg^  j/\AtAu™tu~\uvr\ir\ 

•*^^gr^gr^gr-*sr4>j/grg>fgr.f-*4>^^g*grg»g»gr>J'g>-grgr  t  *g>  •*  ~*>r>**t>t  *t>t*t*T>t***'*t*t't***>***-t  *t  -t^t-t  t-r-t^t^tfj-^**^ 


uj  NOino  omcj-iA  — «>o-^ir»nD— •  g»aorgr-  r or«-gr— «nomrfv\i— •(r-a^cor-gjiAgr ^-^cj*aD-Ovr*^oa»<>f^p*QDC>'^*nr*coo-^fn'nr^o,N*rninr-'(^ivjr^ojf^ 

O  r-ugjouMAgrgr  -J/ ^nt\-<^ o (j*^coa3T3 r-r*r-  ogogD^-u^^u^Au^u^iAg^ina^^grgrg-^g-rnrnrnmrng^gr^gr^^ 

«X  <^fn^rnr»i^rnr^r^i^rO^(^VMa>rgr\jt\if\jrN|rg^rgrsjr\jrsjrgrsjrvjrgrsj  '\jr\jr\jcgpg'"\jf\jr\jr\j  rg  '\jcgrgogrgr\jrgcgrgfv'\|iNjr\jr\j  Ajrgr\jr\jr\j"grgr\jrgrgr\j<-\jrgr'g 

X.  O a*  cor-  ■otn^m'Nj-^ooaiN^in  v^rnrgr-^ocT^cor-oirv^r^rg— <Oi>cop»-  omvtmfNJ^o^^s  OiAvymrg-*oo>cor~  oiAgrrorg^oo^cor-  -oiAgf.-nfNj 

^,  ocA(j*lt  c»  <j*  ^^<j*  o^o^co  co  rococo  coco  .o  t)  a*  r-  r-  p*  I*-  r-  h-  Is*  r-  r-  r*  o  -o  gj  gj  gj  gjNogjgjgjiAiAiA»AiAvAiAtAUMA^r  gr-rgr-j/grg-grgr  •rmiofonmminrt) 

O  ••••••99 •     * , »••••••• 

^oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

o  m 

o  -> 

a  a. 

O  0^g*^Nrxjr^o<\j^^Nr^r-CT'Ou^gfa3N^C^>J'Ono^^g?'co<N^ 

vi  gprnrNjgj^coo^c^rNjg^gor^cornrn^^cDr^f^og^cDgrN^                                               i^coNMr^-^ocof^Na-^ooin  i-+t\jirt  gj  cor»-rsj  -rtAr-rg  -ocnjo"* 

>  — «r-mgr  -oo<- j o^  o  go—nr*  Amcooo^A  g'-'Ag")^— <  or—rgj-  r-  og^r-tArgrno^r-  ruor*-r^iAOao<~gij,tO'"0''Org  rONin«t  ^■•n— 'on  Ti^gr-'joiAcorj  r-o 

UJ  rNJcOU^OrNj^st^^^O'^^O^O^N  Ornr^^^r^orvj^ooor—  ^                                                                                             gr-  T  /^T'lnr^mtn^fiaO  -O'n-^^I/'O  O 

CJ  r— r— cj*<"\j.nLr"^co-o<"og"<j^m-*o  ^^^NON,inr^-*otTiii(r'Q-«rn  ■utJ,,orNJ-Tgj<r'c\jiAro~*  ^r— ^•Acoror».-v.rr--rvgr--rnt7*iA--«r»-  gr  — 4-o-j-— *  r> r- mm  gr 

A.  r»*nm>^o,-<jr-u^g'ptvj  — «^rujr*uA  ^n-ojt^  uu^ni^o^^'w^  omu^^ri^j^^oLPi^ad^r*  ■uu^Jigrg>^«M»\jrg^^^ooo^LrtJsD<u',j-.Jjuj 
grgr^grgrrMr^Pi'inimtgrNj^cgrgrgrNjrNa.^^^^^^^'^ooooooOOOO 


>        • • t 

UJ  OrNiNrnjOO^O-CAO  0'NJC0grOiA~HAOgDg>-g>grr*')rAgrg>-iA  O^O^r^Oebi/>c^^^cO^^-^(7'f^(^^OC*cO«*IO"^<OrOh-Of*>»00,*0-H(\|fOi-4rvJiAr*"0^ 

O  rgoo^r— gDiArnrg^g^aogjiAgrrg— «<j^aogDiAgrntf\j  — <OCn'Dr->o.ym^ui^^'  gr  >r  mrAmmrvjrgrg  — tr-i^uO^-^coaDr-r-  gjiAgr  T<*1CV- <  ^ -^ —« ,-n -g  c  g; 

i:  OJ'CON  OLri-TrorNi— *oc7,icn^--OiAst^rNJ'^0<7s  x»N-gDa^  ^fmrsj-to^coN  OiA  j-mtNi^HOO^cor*-  gpiA>jT^fg^O<7,cor--gjtA  g*  «vnj  — <  o '7*  co  r—  Oia -t  f*"!^ 

<i  o^o>^O^O^v*C7*t^crcT'^*oococooacoa3a3co^r— r^r^r— r-r^                                           -OsO  OAiAiA»AiAlAtA,A1AiAg>g,,grgr  4-  gr  g^  gr  sj--jTO«irnmrr\mrt^  n 

O       ••••••••■■ .............. •••••■• 

o  -< 

O.  Of 

ft  ftft«««ftftft«ftft<ftft«ft«««&ftftl»ft^«««ft«ftftft»ft»ftftft^ft«««ftft«-«ft«ftftftft*ftft«ftftft«««ft«ft«« 

cj>  rgcowOcouco^grgrg3g}*x>rgoAj<Ngrgro^cog>xa}gr^jT>ogJ^                                                       sro^J,ojrM-oor«.ococo*T  sOJ3rNi^>r-t-tOO^^OM-OJ" 

%/>  f~grr\jr\j<-»ao^v~>o<7<r*coiAr-iA^g3gr  colt-co— «gj  g3vr<^o^r«a^'n'X>co~-<cocorgr-oi_>  j-  ^u^cNjgrrgaao^-*— «u-\'<ig'r\jsugrrrvM— « u~v  <-g  gj  cm  o  o **\j  u^ r-  r*  -o 

^  rvrvjr^^^^^t/tu^^4Drno^iJ>u~tg^-4JfN*Of^co'NU~\--^^-«u-\^-<iJO  u^Jfg  j^u^i^ot^orgcoocoo^J^'^f^^w^oO^J^^f^nrALrLO-r'g^^a^nirNjxjsj-^ni^ 

uj  ONUiNO^N^*nf^^f^rnNO^,w^ij'Qinin^ini^>ocoN'00^0'^ocuo^r-in>rvr>orn'>j<Hr^  ^roocr  r»->. oj so <j y^  iri.-gcy-^^o'-o-o'Ni'-A'-i-'or—  j«i 

Q  /^r^ro<^^rsjr^rsjcoNfOc04DNOO^coa^rsjr\jojrnaDOiA-^a)in^rNjrnrO'^^^                                      r\j,o— 4Ln^i\i^>'0  j-aorgr^-ng  o-^n^  o^- trMoxOi-g 

f^  rO  n^ -^  r^  r<^  rO  ro  ro  r^  rO  n^  ^  fO  ro  rn  rn  rn  .^  rt^  rr,  rr>  ro  <^  ^  i^ -^  f^ 


> 


Q     ^Hrgrs;rrK-^r*irnOf,grvJ— <OCJ^cof^r»-g3^-0  -oltmalTi  fgf^tin  0'Or*7«OvOOO^^-<'^-<,vioj(\jrg^Hr-4^oo"N  ogf\J^oco  r*r*--o>o  o  ^o-o^O  -OiA 

iT  o<7k<B^OingffOfNi^^OO*oo^g3irt<f,^c^'-^ochoofs-g3U"»<fcorg]i-<oo>eor*  >omgf  <noj-HOO*uoN^irkg|Tnco>^oovcor,*gDu^gh<,n<v— ^o^oon-  -om^rotN 
<S  O^O^^O^Q^C>^(7vC7N^cocooOCOUDa3cocococOfs-'s-r-r*r,»N^f,,"r,-r^  ^^g^gOg^OOg^gDOiAiAiAiAiAiniAiAtAi^-rg^Nj'  gv-^,g^g>gr^grr^r^^r^r'>mf<lro 
o      •••••••••« t ••••••••••••■■••■••••■••••••••••••■■•«•••••■■ 

^oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


•X 

<D«-of\oooooo^j"OOOTj-oi'jx-Tr\j>r -ooj-cosr^o  -rrsio  uoo-u'o-O'M^W'j'wro^o-oN  oaioajaj^r-ro^j  o-u-ro-o^-r-r-rorgo 
Oci'Miru  ru^^o^  mo'nI-t  -nt^_0'\iir\<x>ijM\jujo  u  1*1  o  r-<rgiJ,'O-t"^0u'»tr<^iO  *jco<\j  yy-r  -«vnuvrv  r  aoro^fg  «t:mnu»j- rvjnjajgjg-aj-TjT./* 
0^^iAONf\ITi^i>iJvO  U'7» -T  TJ"'0<NJ   *uujr'|inO's'iJ.'^'-<-*'-<  OOON^O^O'JJOaD-O'MO^JNin  0-*?0-0«f  U  <X)  M  ■  A  O  ^  —<  r—  DfNJ- «iAC       J  O 

(j^rNjiNjg<rx)N>r'Trrior--iju»o-rjon^srgi-<"ur^-o--*-0'XJ^g"  juro^o^  o^HNiMnotn  r— <'/>-r— «-- i-ro-^no^  oia— •-o^ir-r-r'v^'-ro)  ^,,  j— «  -\j 

r*^0  g"  >x>~*  >jj—*r~i*)^j  -o-n— *co•4Ju^  ^m<^o,-T'VJ-'r,,--ugJ-OlJ'>lJ^  oiaum/^ia  ofs»co^-<g*'fs,»-H  •y'JO'NJQ,-U'**|n<"»^«r'DOf,HiO'rtN*-<-T'**  a^r^o  o^r- 
y^^cof^r^gjxJ^j^a^^vrgfrn^ifnrnrornrn^rgiNiryrgrgrvjrsjcNjr^  -J/iAia  <j^*f^co^O"-*^gtNJtvrrmni^--r  g*'AtA  J^u"MAtA  sr  j-  rn 

rggr^Jfvj'Njr>jrsjr\jr^i%jrg^r\jr»jrNj^rgrNjrgogfvfsj'Ni,^rg'^'^'Nf^^'^^ 


Uj  ^lr*-HJ^frlrr^^H^^i^r-g>C^'^  JOO--«  O  iA  O  — « ^  lA  CO  CO  r- g."  COCrCT*  OOO— *r^-^r\j  aOcoaOOlAOiA  >CD^J*0^~*-*■>'"Og,•  4->J"X)^J,gJCOr\i-OcOOCO^^t7'Or-«CO 

c_)  ovAojaogj^gr"  <^(>jrNjrg^^^a3ajcDa5aDP^f^a3a3CO-i3CC^coflX)JN\JNC^^O^cggr>0                            j>  jd—*  j-r^<.-rNj  gj- gjr*-a30-— ^igr -A or^-OiAfA rg^in 

<t  iAgrgr<^^cnro<^^rnrOror^pg<NJ^gi\irsjrgrg^jrsjr\j^pgrgrNjrg<^                                                               ^irnnmo-ocor-r-f-f-^^i         ~ococcro''-ococor»r» 

S:  CO  U3  CD  CO  O?  CO  CO  tO  a>  03  CU  T>  CO  03  CO  03  CO  ^  <X)  ^  U^  (^  CO  T3  CO  CO  CO  T)  GO   O^J  0Q00CO  OJCOCOOJCOCOCOajCUCOCOajaOajTJ^  I^COWCOCOCy  CO"A»^      J»-CU              '^O  00  co  co  co 

3L  oo^u^Ngju^g^f^rNj^ocrajf^g^u^Nrmr^^Oc^cup^gDU^grnrNj^o^^f^ 

^  u^y^^ijN^vJN^o^^^ajajcoajcococoJj^ajr^r-r^r^r-i^r^f^r^                                                                                                                     -  <r  jncnmnminnni 

o  »«#«§»•»•■•»#•••••■•••■•••••••••••••§•••••■••••••»■•••••••■•••••••*•• 

hoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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.  .  COMPUTER  OUTPUT  FIVE -THREE 

Aggregate   Mean  Absolute    Deviation    (in    cubic    feet) 

OOOOOOOOOOCJOOOOOO 

iuuiujiij*ujtuujiijiijujuiuiiuuiiuuiiu 

•o  o           *           ooioo^co(M<oo>rcor\j>oo 

0  ♦  •               •               •               •               t                •               «                •               »               •                t                •               •               •               t               •     o 

tu»»««»+»*»-»»«t)l»»*»»»*»»«»0»*4.»#»*»»«*«»#*»«»+»»»»*»»+»»»»»««»««««... 

«->♦  ♦                                                                                                                                                                                                                                                                                                      *o 

o«  ♦                                                                                                                                                                                                    »o 

••  ♦                                                                                                                                                                                •  . 

♦  ♦  • 

♦  .  *                                                                                                                                                                                ♦ 
-l»  ♦  ♦                                                                                                                                                                                                »-• 
o»  ♦                                                                                                                                                                                      «o 

1  #  ♦                                                                                                                                                                             •  I 

ui*  *                                                                                                                                                                                                                                                    a  ,jj 

o*  ♦♦                                                                                                                                                                                             *o 

o»  ♦                                                                                                                                                                                             *o 

•  »  ♦                                                                                                                                                                                   *  • 

#4  • 

♦  ♦  ♦ 

-«»  ♦                                                                                                       .                                                                                          »  -I 

©•  ♦                                                                                                                                                                             »o 

1*44  «  I 

l"»  4                                '                                                                                                                                                                                                      « tu 

o*  ♦♦                                                                                                                                                                                 *o 

o»  ♦                                                                                                                                                                                           »o 

••4  «    • 

•o«  *                                                                                                                                                                              » o 

♦  ♦  .              » 

♦  4  ♦ 

-l»  ♦                                                                                                                                                                                                                                                                                                                  •  -I 

o»  ♦                                                                                                                                                     ,  __                »o 

I  •  ♦                                                                                                                                                                 •  I 

UJ«  »                                                                                                                                                                                                                                                                                                             IHJJ 

o*  ♦                                                                                                                                                                                            *o 

o«  ♦                                                                                                                                                                       »o 

•  »  4                                                                                                                                                                                                                       *    • 

■*■*♦♦  *  •»• 

«•  ♦                                                                                                                                                                              •» 

♦  ♦  ■»• 
-<«  ♦  ♦                                                                                                                                                                          *^< 
o*  ♦♦                                                                                                                                                     *o 
I  *  ♦                                                                                                                                              *  I 

OI»  ♦                                                                                                                                                                                                                      *w 

o*  ♦                                                                                                                                                                +  o 

O*  4                                                                                                                                                                                                                              *0 

•*  ♦♦                                                                                                                                            *  • 

ev*  ++                                                                                                                                                          *w 

♦  ♦♦  * 

♦  ♦  ♦  + 

00»  r                                                       ♦                                                                                                                                                                                                    *<D 

o»  ♦                                                                                                                             »o 

I  »  ♦                                                                                                                                            ■»  < 

0  +     .••.«•••••••••«••«••••«•••♦     •  +     ......................     I....     ..t.....+   +   <3 

</<»  ♦                                                                                                                                » >^ 

.*  ♦                                                                                                                                                                  *r    • 

in*  ♦•  ♦                                                                                                                                  #m 

1  *  ♦♦                                                                                                                                             *  I 

♦  ♦♦  ♦ 

o«  ♦                                                                                                                 *o 

I  *  ♦                                                                                             ♦  I 

ID}  ♦                                                                                                                                                                                                                                                            »UJ 

o+  ♦                                                                                                               *o 

o*  ♦  ♦                                                                                                                    «o 

•#  ♦                                                                                                       *  • 

rg»  ♦                                                                                                                                                        »<M    /— n 

l»  ♦♦                                                                                                                                            *  I       nj 

-l»  ♦                                                                                                                                                                                                                                                                                                         *-•          fc 

o*  ♦  ♦                                                                                                                                                *o     £ 

I  *  ♦       ♦                                                                                                                                                                                                                     *    '        rrt 

o*  ♦     ♦                                                                                                                                                                           +°     W) 

o*  ♦                   ♦                                                                                                                                                                                               *°   yj 

•*  .       4                  4                                                                                                                                                                                                       »    • 

>r*  ♦    ♦                                                                                                                  *■? 

I*  ♦+                                                                                                            *  i    ji 

-.«  ♦                                                                                                                           *-    £ 

oit  4  4                                                                                                                                                           » o      _, 

uj«  4*                                                                                                       .                                                                                                 *^     +-> 

s: 

.*  .                                    ♦  ♦ 

?:  *                                                       »t  o 

♦  ♦                                                        ♦   u 

Ofl  ♦ 

I   o  * 

urn  * 

°*  I 

O*  4 

oo»  * 

■  »  + 


»0       V) 

*   .     d 


»  •  bO 

it  Ol  ^ 

♦  O  »-i 

«  o  .H 


;Q  o 

»o    O 
♦  o  CO 


o*  * 

o»  * 

ill*  * 

S»******#»**«»*»*»***«**»*i-»***»**+***','*********#*+*****************0 

-<oooo<->ooooooooo<jo 

UJa(UjUjujaiuJuJ"-'iuiU"J  ujaiaiuj 

_,Oi>  oirvtnrgocor~in^  fsioi^r- 
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COMPUTER   OUTPUT    FIVE -FOUR 


— «/\  g*.-«  J"  — «-<rvj   f  0'JJ"0-<^n.">  OO  JT  Jv^-<0  ^O^l^*"'O^OOuJfjyi  J'^TO!'  It—  OX  C7""gOiA  -O  N"".^©  J"^   0^"^*^^- «-7>r»*  t)^^ 

p*<Dao^OAj'>gg*iO'--o-*on3*-^-^.r-r,-r-r-r*-<oo  — ot)^  o-oia'AiAia  j-vini/vAtAu-MAiAtAiAiAiAiAiA-r  g*  ^-r  j*g>g>g>g>g*g>g>g>g>>rg>g>g>g> 
\r\*f\%r\m -o  *o  o^r»«f~^^ir»ir*ir..n  j"\irw\u"^m»/\u"\u^.o  o*Ag"*«*g~»iA'AiAg-»iA<AiAj"\iAtAu"w*g-\  A^^g^ig^u^u>tf^g^unj>u*vA'AiA«AiAirw'»*AiAiAvni* 


-o^'T^'O^-f  -TKsoo>f  -*r\j.>m----mf\j^«  o^  TtN»-r-o'fc-r-r^^r^^r<3r«-r«--n— «^^oor^j'MuriA)'*irf— r"-'*^'Min— «oxia<"*">Oia  -riAo-*orriu-» 

•— •j,<iy*lT'N*'T3UVJji^aj>T0D»J'_Jnj(j>ri^jr*^v/>  xj  — on^O'M-O  —  ^uiaj  <j— <o— *^roaj^ja\jj  ^^^r-rni/-  ^»J«u\*M(/>*j-r— «O0  — «*M  T  ™-  Jf^Tj^y 
aj  iNff*  g>  -TiNJ— «CO  gf  O-^i^JC^'-jr^-iMrvi—i^H— t— imr-mr\j^-*LT:n^  4* -*  O*  CU  .J  CO  CO  JO  V  0s  O—"^''  •nni'NJM^iHQOO^  <J  >/>C7,U'  LrU"1  '.f  J*  l^OO-^'-'W 
pgi'tVMJ'in  ©•0-0330^-u"\«~ir'>r\j.\)  >j<\j  Njrvj-vj'-grNjf^  >r  ,j  1~grg~«—t.—ioOOOOOOO—*— *—♦—•—«—«'-*--«—«—«—«"'— 'OOOOO  OOOOQ-*-^— <— •-* 


OOOOOOOOC  OOOOO  OO  OOOOOOO  OO  OO  OO  OO  OOOOOO  OO  OOOOOO*  JO  OO  OOO  OOOOOOO  OOOOOOO-* 

I  I  I  I  I   I  I  I   I  I  I  I   I  I   I  I   I  I  I   I  I  I  I  I   I  I   I  I   I  I   I   I   I  I   I  t  I  I  I  I   I  I   I   I   I  I   I   t  I  I   t  I   I   I  I  I   I   I  I  t   I  I   I  I  1   I  I 
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